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SUMMARY 

The  development  of  a mathematical  model  for  fatigue  crack  growth 
under  arbitrarily  varying  cyclic  loads  was  initiated  under  work  previously 
supported  by  NADC,  as  documented  in  report  NADC-76076-30.  This  model  utilized 
the  inclined  strip  yield  superdislocation  representation  of  crack-tip  plasticity. 
While  in  good  agreement  with  experimental  results  in  predicting  steady-state 
crack-growth  rates  for  constant  AK  loading,  the  model  was  found  to  be  Inadequate 
in  other  aspects.  Specifically,  the  mean  stress  (R  effect)  was  not  properly 
reproduced.  Also,  because  the  model  retained  one  superdlslocatlon  pair  for  each 
load  cycle,  the  computational  procedure  was  quite  cumbersome. 

Steps  to  overcome  these  deficiencies  are  described  in  this  report. 

These  were  undertaken  using  a simpler  model  in  which  one  superdlslocatlon  pair 
was  used  to  represent  the  residual  plasticity  left  behind  by  the  growing  crack. 

In  addition  to  providing  results  for  constant  crack  growth  rates  that  compare 
well  with  experimental  results,  this  work  has  given  some  basic  insights  into  the 
problem.  It  has  also  indicated  how  improvements  in  the  model  formulation  for 
general  load  sequences  can  be  made.  Tests  of  the  model  show  that  the  main 
features  of  the  experimental  observations — crack  closure  Included — can  now  be 
reproduced.  While  not  yet  to  the  point  where  crack  growth  calculations  for 
engineering  purposes  can  be  routinely  performed,  it  is  believed  that  the  model 
development  is  well  on  its  way  to  this  goal. 
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PREFACE 

This  research  program  has  been  conducted  by  the  Applied  Solid  Mechanics 
Section,  Structures  and  Mechanics  Research  Department,  Battelle's  Columbus 
Laboratories,  Columbus,  Ohio.  The  work  was  performed  under  Contract  No.  N62269- 
76-C-0316.  The  program  was  administered  by  the  Air  Vehicle  Technology  Department, 
Naval  Air  Development  Center,  Warminster,  Pennsylvania.  Mr.  Paul  Kozel  provided 
the  technical  liaison.  This  report  summarizes  the  work  performed  during  the 
period  from  July  1,  1976  to  June  30,  1977. 
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1.  INTRODUCTION 


A vast  amount  of  experimental  and  empirical  analysis  work  has  been 
done  on  fatigue  crack  growth  over  the  years.  The  field  Is  sufflcently  well- 
documented  that  comprehensive  accounts  are  now  available  In  basic  text  books — 
for  examples,  see  Broek  [1],  Hertzberg  [2],  and  Rolfe  and  Barsom  [3] — and  In 
various  popularized  accounts  [4-6].  Yet,  at  present.  It  Is  nevertheless  true 
that  the  fatigue  process  Is  simply  not  well  understood.  As  a consequence, 
there  are  many  situations  In  which  fatigue  crack  growth  behavior  cannot  be 
quantitatively  anticipated.  Perhaps  the  most  outstanding  example  Is  the 
current  Inability  to  predict  crack  growth  In  cyclic  loading  sequences  that 
Include  load  level  excursions.  To  treat  this  and  related  situations  properly, 
there  Is  a distinct  need  for  a quantitative  predictive  method  arising  from  a 
fundamental  mechanics  point  of  view.  This  report  describes  a model  which  Is 
aimed  at  fulfilling  this  need. 

A predictive  model  for  fatigue  crack  growth  under  variable  amplitude 
loading  Is  faced  with  two  conflicting  requirements.  Thf  first  Is  that  the  crack 
tip  deformation  must  be  accurately  predicted  In  the  presence  of  the  residual 
crack  tip  plasticity  generated  In  previous  load  cycles.  This  Is  necessary  for 
the  model  to  be  realistic.  The  second  requirement  Is  that  the  model  must  be 
simple  enough  to  permit  a great  many  load  cycles  to  be  calculated.  This  Is 
necessary  for  the  model  to  be  useful  for  practical  applications. 

A mathematical  model  that  Is  an  effective  compromise  between  these 
requirements  was  Initiated  In  previous  work  supported  by  NADC  [7-9].  The  model 
Is  based  upon  the  Inclined  strip  yield  superdlslocatlon  representation  of  crack 
tip  plasticity.  This  gives  rise  to  a ntimber  of  extremely  desirable  features: 

(1)  It  admits  load  levels  that  can  vary  arbitrarily  from  cycle  to  cycle,  (2) 

It  automatically  takes  account  of  the  crack  closure  phenomenon,  and  (3)  It  Is 
capable  of  predicting  crack  growth  Increments  solely  from  readily  available 
mechanical  property  data.  Thus,  of  most  significance,  no  disposable  parameters 
connected  with  specific  load  spectra  are  required  to  perform  crack  growth 
computations  using  the  model. 
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Predictions  made  using  the  Initial  version  of  the  model  were  not 
completely  satisfactory.  This  work  was  undertaken  to  remove  the  deficiencies. 

To  test  the  results,  specific  attention  was  placed  on  the  steady  state  crack 
growth  rates  resulting  from  simple  constant  amplitude  loading.  As  shown  In 
this  report,  a closed  form  relation  can  be  extracted  from  the  model  for  these 
conditions.  It  Is  found  that  the  result  closrely  resembles  the  semlemprlcal 
fatigue  crack  growth  relation  developed  by  McEvily,  et  al,  [10-11].  McEvlly's 
relation  has  been  found  to  be  In  quite  good  accord  with  experimental  data  for 
many  different  materials  under  a range  of  loading  and  environmental  conditions. 
Consequently,  the  correspondance  is  taken  as  confirmation  cnat  the  model  Is 
now  in  good  agreement  with  well  established  experimental  results.  Other 
experimental  comparisons  are  provided  which  reinforce  this  conclusion. 

The  agreement  between  the  crack  growth  rate  predictions  of  the 
present  model  formulation  and  constant  amplitude  fatigue  crack  growth  data 
constitutes  the  most  significant  result  contained  In  this  report.  Nevertheless, 
it  is  Important  only  Insofar  that  it  demonstrates  the  correctness  of  the  basic 
approach  taken  In  this  work.  It  Is  but  one  milepost  on  the  road  towards  the 
development  of  a general  predictive  technique  for  fatigue  crack  propagation 
under  variable  amplitude  loading. 

The  work  presented  in  this  report  demonstrates  that  the  model  Is 
fundamentally  sound  and  strongly  suggests  that  It  Is  capable  of  providing  an 
engineering  tool  suitable  for  fatigue  crack  growth  calculations  In  service 
conditions.  However,  at  the  present  stage  of  development,  the  model  Is  too 
cumbersome  to  be  used  routinely  for  complex  load  spectra.  Further  work  Is 
required  to  develop  suitable  representations  of  the  residual  plasticity 
associated  with  arbitrary  loading  histories.  The  simple  model  developed  In 
this  report  for  a uniform  cyclic  loading  shows  how  this  can  be  accomplished. 
This,  together  with  an  approach  that  will  Incorporate  environmental  effects 
into  the  model,  will  provide  a broadly  applicable  fatigue  crack  growth 
predictive  model. 
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2.  PRELIMINARY  DISCUSSION 

Fatigue  crack  growth  predictions  made  using  the  preliminary  version 
of  the  model  were  found  to  be  at  variance  with  certain  well  established  experi- 
mental trends.  In  subsequent  work  undertaken  to  remove  these  deficiencies, 
specific  attention  was  placed  on  the  crack  growth  rates  resulting  from  con- 
stant amplitude  loading.  A large  amount  of  experimental  data  exists  for  this 
case  that  can  be  used  for  comparison  with  the  model.  Also,  a closed  form  crack 
growth  relation  can  be  extracted  from  the  model  to  facilitate  these  comparisons. 

In  fact,  when  this  Is  done.  It  Is  foimd  that  the  result  approximates  the  fatigue 
crack  growth  relation  developed  and  verified  by  McEvlly  and  his  associates  [10,11]. 

In  this  section,  McEvlly 's  semlemplrlcal  crack  growth  equation  is  first 
given.  Following  this,  the  basic  equations  for  the  model  are  described  together 
with  the  procedures  for  obtaining  relations  appropriate  for  linearized  small- 
scale  yielding  conditions.  These  discussions  are  Intended  to  provide  background 
for  the  material  presented  in  subsequent  sections  of  this  report. 


McEvlly 's  Semlemplrlcal  Crack  Growth  Relation 


McEvlly  and  his  associates  have  examined  the  various  power  law  and 
other  relationships  that  have  been  advanced  for  constant  amplitude  cyclic 
loading  In  the  light  of  experimental  data  for  a wide  range  of  materials,  load 
levels,  and  environments.  From  this  work,  and  also  partly  guided  by  theoretical 
considerations,  they  have  evolved  a new  crack  growth  law  of  great  generality. 
This  law  can  be  written 


dH 


AK  - (AK) 


T 


1 + 


(1) 


where  E,  Y,  AK,  fL,  , K , K . , and  R have  their  usual  meanings,  here  and  else- 
ic  max  min 

where  In  this  report,  as  follows: 

K • stress  Intensity  factor  at  maximum  load 
max  ' 


^mln  ~ Bbress  Intensity  factor  at  minimum  load 
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E elastic  modulus 

Y » the  tensile  yield  stress 

» the  plane  strain  fracture  toughness 

AK  » K - K . 

max  min 

Also,  A Is  an  empirically  determined  environment  sensitive  material  parameter 
and  (AK)^  Is  the  threshold  value  of  AK  below  which  no  crack  growth  Is  possible. 

Values  of  the  coefficients  In  Equation  (1)  for  several  different 
materials  and  environmental  conditioi “ are  given  In  Table  1.  A comparison 
between  the  predictions  of  Equation  (1)  and  experimental  data  on  T1-6AJI.-4V  at 
two  different  R values.  In  air  and  In  vacuum,  which  exeiiq>llfles  the  accuracy 
of  this  relation.  Is  shown  In  Figure  1. 

Equation  (1)  contains  a stronger  R dependence  than  It  might  appear 
since,  according  to  McEvlly,  et  al. 


(AK)^  - 


(2) 


where  (AK)^  is  the  threshold  value  at  R “ 0.  The  equation  also  exhibits  a 
strong  yield  stress  dependence.  However,  the  parameter  A can  depend  on  yield 
stress.  Specifically,  In  Inert  atmospheres,  A Is  proportional  to  the  ratio 
Y/E.  Inserting  this  Into  the  above  equation  and  rearranging  then  gives  a 
result  of  the  form 


While  this  evidently  Is  Independent  of  yield  stress.  It  should  be  kept  in  mind 
that  It  Is  not  a general  result,  but  holds  only  for  tests  conducted  In  a vacuum. 

The  most  significant  point  about  McEvlly's  formula  Is  the  fact  that. 

In  contrast  to  the  conventional  power  law  formulation,  the  AK  dependence  Is 
quadratic.  This  Is  Important  because  the  mathematical  model  being  developed 
here  Is  associated  with  a crack  opening  displacement,  strain  energy  release  rate 
point  of  view  and,  as  Is  well  known  in  fracture  mechanics,  these  two  quantities 
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TABLE  1.  VALUES  OF  COEFFICIENTS  IN  SEMI- EMPIRICAL  FATIGUE 
CRACK  GROWTH  RELATION  OF  McEVILY,  ET  AL. 


Y 

E 

Kc 

ksl  In** 

A 

(dimension- 

less) 

(AK)o 

Material 

(ksl) 

(ksl) 

ksl  In** 

Ti-6A1-4V 

Forging  In  Air 

120 

Titanium  Alloys 

17x10^ 

65 

.05 

6 

T1-6A1-4V 

Forging  In  Vacuum 

120 

17x10^ 

80 

.02 

12 

T1-6A1-4V 

Sheet  In  Air 

140 

17x10^ 

100 

.08 

5 

2024-T3 

57 

Alualnum  Alloys 

10.8x10^ 

100 

.03 

5 

L64 

52 

10.6x10^ 

60 

.18 

6 

RR58 

60 

10.7x10^ 

21-24 

.04 

4 

A514-B 

129 

Steels 

30x10^ 

200 

.02 

8 

X65 

65 

30x10^ 

150 

.02 

5 

En30a 

Uneobrlttled 

107 

30x10^ 

200 

.026 

7 

En30a 

Embrittled 

107 

30x10^ 

50 

.026 

7 

316L 

41 

Stainless  Steel 

30x10® 

90 

.01 

8 

I 


5 


) 

i 
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FIGURE  1.  COMPARISON  OF  EXPERIMENTAL  CRACK  PROPAGATION  RATES  FOR  Ti-6AJt-AV 
FORGING  IN  AIR  AND  VACUUM  FOR  TWO  DIFFERENT  R VALUES  WITH 
PREDICTION  OF  EQUATION  (1);  DATA  OF  McEVILY,  ET  AL  [10-11] 
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are  related  to  the  second  power  of  the  stress  Intensity  factor.  Thus,  as 
shown  In  the  following,  the  model  will  produce  a result  similar  to  that  given 
by  Equation  (1).  One  difference  is  that  the  dependence — which  may  be  due 

to  a crack  growth  mechanism  occurring  at  high  load  levels  that  differs  from  the 
mechanism  that  predominates  in  the  near  threshold  to  moderate  load  regions — is 
not  taken  into  account. 

Basic  Equations  of  the  Model 

A coi)q>lete  derivation  of  the  equations  governing  the  superdislocation 
representation  of  crack-tip  plasticity  have  been  given  in  references  [7-9].  For 
present  purposes  it  will  suffice  to  recapitulate  the  key  results.  First,  for  a 
tensile  loading  applied  in  the  direction  normal  to  the  crack  plane,  for  plane- 
strain  conditions,  the  shear  stress  t acting  on  each  superdislocation  is 


+ : r 2 ) , n = I , 2 . ...  M , (4) 

''  8ir(l-v  ) j=.i  J J"  J” 


while  the  stress-intensity  factor  K is  given  by 

''  ■ I”  - efe  ‘I'il  • 

where 

o ■■  applied  tensile  stress 
a “ half  crack  length 
M ~ number  of  load  cycles  imposed 
bj  ■ strength  of  the  dislocation  created  in  the  j—  load  cycle 


with  the  remaining  undefined  quantities  in  Equations  (4)  and  (5)  being  known 
functions  of  the  dislocation  positions.  These  are  given  in  terms  of  the  conplex 
variable  representation 


'j  'J 


19 


1 . 2 


M 
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where 

e ~ angle  between  the  crack  plane  and  slip  plane  (constant) 

a.  = crack  length  at  the  time  that  the  j—  dislocation  was  emitted 

th 

■ distance  from  the  crack  plane  along  the  slip  plane  to  the  j — 
dislocation. 

It  Is  considered  that  the  material  provides  an  Intrinsic  friction 
stress  which  resists  the  movement  of  the  superdlslocatlon.  This  parameter  Is 
denoted  (for  Internal)  and  Is  related  to  Y,  the  tensile  yield  stress  by 

■ eY  , (6) 

where  e Is  a numerical  constant.  The  equations  of  force  equilibrium  that  are 
the  basis  for  the  mathematical  model  are  formed  by  setting  the  right-hand  side 
of  Equation  (A)  equal  to  for  each  of  the  M superdlslocatlons  In  the  problem. 

An  e]q>resslon  for  the  crack-opening  displacements  can  also  be 
determined.  Again  omitting  the  details,  the  result  Is  found  to  be 

,2  , M 

V - Ao  (c®  - x®)^  + sin  0 E b d , (7) 

E j.l  J J 

where  v ~ v(x)  Is  the  crack-opening  displacement  at  a position  on  the  crack  face 
X ^ a.  In  Equation  (7) , x Is  measured  from  the  center  of  the  crack  and  d^  Is 
again  a function  of  the  superdlslocatlon  positions. 

For  most  fatigue  crack-propagation  problems,  a linearized  or  "small- 
scale  yielding"  approximation  Is  valid.  For  a plastic  zone  that  Is  represented 
by  a single  superdlslocatlon  pair,  the  governing  equations  are  for  the  force 
equlllbrlm  of  the  dislocations  (only  one  needs  to  be  considered  because  of  symmetry) 
and  the  singularity  canceling  equation.  The  unknowns  are  the  superdlslocatlon 
strength,  b,  and  the  position,  1.  Setting  t ■ In  Equation  (A)  and  using 
Equation  (6),  It  is  found  that 

i 
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eY  - ah  (4)  + [g  (A)  + k (A)]  , 


and 


^ 

8ir(l-v  ) 

In  small-scale  yielding  where  1 « a,  the  functions  appearing  In  Equations  (8)  and 
(9)  are  simply 


f(f)  =6  sin  e cos  I 

g(-t)  “ - -J  (1  - cos  2 e) 

1 ^22N^  6 

h(f)  = ^ \~l)  ® '2  ® ® 

k(i)  “ " ^ cos  20  . 


(10) 

(11) 

(12) 

(13) 


Substituting  Equations  (10-13)  Into  Equations  (8)  and  (9)  and  eliminating  b 
gives  the  following  relation  for  the  superdislocation  position 


I 

a 


! , . e 1 i 

a 

1 (‘ 

1 1 

a I 

\ -j  s *ri  c cos  A p.  1 

( sin  0 cos  ' 

1 

() 

fi-T 

Y sin  0 coo  0j 

-2 


(14) 


In  this  work)  the  particular  angle  6 selected  was  that  which  maximizes  the  extent 
of  yielding  (equivalently » the  plane  of  maxisnim  shear).  From  Equation  (14),  this 
la  e - cos"^(l/3)  - 70.53”. 
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The  superdislocation  will  not  give  precise  Information  about  the  size 
and  shape  of  the  plastic  zone.  However,  this  Is  unimportant  because  It  accurately 
represents  the  plastic  deformation  occurring  at  the  crack  tip.  Consider  the 
prediction  of  6,  the  crack-tlp-openlng  displacement.  This  can  be  written  as 

6 - 2b  sin  e . (15) 

Upon  solving  for  b and  substituting  Into  the  above  for  plane  strain  under  fixed 
load  for  small-scale  yielding  conditions.  It  Is  found  that 


3 cos 


1-  V 


0/2 


sin  0 cos  Y - 


sin  0 cos 


EY 


1 - ® cos  0 


(16) 


where,  for  an  Infinite  medium,  K - a / utt  . Hence,  6 - a — , where  o - 0.A33  for 

Ell 

e “ If -/T aad  v • 0.3.  A comparison  with  the  results  obtained  from  a variety  of 
both  approximate  and  highly  rigorous  calculations  of  the  crack-opening  displace- 
ment are  given  In  References  [7-9].  These  results  verify  the  fidelity  of  the 
superdlslocatlon  representation  In  the  area  of  main  concern  here — In  Its  pre- 
diction of  the  crack-tip  crack-opening  displacement.  Further  validation  of 
this  approach  Is  given  In  the  following  section  of  this  report. 
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3.  COMPUTATION  OF  CRACK-OPENING  DISPLACEMENTS 

Two  areas  have  been  Identified  from  which  the  inadequacies  of  the  model 
developed  initially  have  stemmed.  These  are,  (1)  the  treatment  of  a residual 
dislocation  in  the  unloading  phase  of  the  load  cycle,  and  (2)  the  role  of  the 
singularity  canceling  equation  in  determining  an  equilibrium  state.  The  key  to 
improving  the  model  in  regard  to  these  two  areas  has  been  the  development  of  cor- 
rect expressions  for  the  crack-opening  displacements.  The  procedure  and  results 
are  completely  described  in  Appendix  A of  this  report.  For  present  purposes  it 
is  enough  to  say  that,  although  these  results  were  derived  in  the  previous  work, 
because  no  use  was  made  of  them  earlier,  their  deficiencies  were  not  recognized. 

It  is  now  certain,  however,  that  the  expressions  for  the  crack-opening  dis- 
placements, which  now  play  a crucial  role,  are  correct. 

When  a maximum  load  is  Imposed  in  the  model,  the  nascent  superdislo- 
cation's strength  and  position  are  determined  by  specifying  that  (1)  the  superdls- 
locatlons  are  in  force  equilibrium  and  (2)  the  crack-tip  singularity  is  abolished. 
This  is  just  the  same  as  in  previous  work.  However,  when  the  load  is  reduced, 
where  previously  the  strength  was  modified  while  the  position  was  held  fixed, 
now  the  strength  is  fixed  (or  moderately  reduced)  while  the  position  is  Increased. 
The  reason  for  the  change  is  based  upon  the  effective  representation  of  the 
crack-face  displacements  thereby  obtained. 

A series  of  elastic-plastic  finite  element  confutations  were  performed 
to  determine  the  exact  shape  of  the  crack  faces  in  a load-unload  sequence.  A 
complete  description  of  the  finite  element  method,  the  confutations!  procedure, 
and  the  results  obtained  are  given  in  Appendix  B of  this  report.  Typical  results 
are  given  in  Table  2.  The  values  of  the  crack-opening  displacement  at  the  center 
of  the  crack,  for  heuristic  purposes,  are  shown  in  Figure  2 for  several  load- 
unload  sequences.  The  idea  was  to  see  if  this  information  could  be  used  to  test 
(and  perhaps  infrove  ifon)  the  superdislocation  representation  of  the  residual 
plasticity  in  cyclic  loading.  In  fact,  this  procedure  has  revealed  a way  in 
which  this  could  be  done. 

The  set  of  calculations  shown  in  Figure  3 demonstrated  the  necessity 
of  modifying  the  superdislocation  position  during  unloading.  In  these  calcula- 
tions, the  superdlslocatlon  generated  at  a maximum  applied  stress  of  29,780  psl 
(roughly  equal  to  half  the  yield  stress)  was  allowed  to  retain  its  full  strength 
and  was  kept  at  the  same  position  while  the  applied  stress  was  relaxed  to  zero. 
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TABLE  2.  CRACK  OPENING  DISPUCEMENTS  FOR  2024-T3  ALUMINUM  CENTER  CRACKED  PANELS 
FOR  A LOAD* UNLOAD  SEQUENCE  OBTAINED  FROM  AN  ELASTIC-PLASTIC  FINITE 
ELEMENT  ANALYSIS  (PLANE  STRESS) 


■ • 1.18  Inches  W - 6.0  inches 

Creek  Opening  Dlsplsceieent  (Inches)  At  Various 
Data  Set  Applied  Positions  On  Crack  Face  (Inches  From  Center) 

Deslg-  Stress 

nation  (psl)  x ■ 0 x ■ .383  x » .709  x - .774  x ■ 1.063 


<•> 

Maxlnua 

Applied  Stress  • 

29,780  psl 

Loadlnx  Phase 

1 

15.190 

.00322 

.00311 

.00257 

.00231 

.00122 

2 

20,050 

.00433 

.00418 

.00349 

.00315 

.00181 

3 

24,920 

.00551 

.00533 

.00450 

.00408 

.00253 

4 

29,780 

.00679 

.00659  . 

.00562 

.00515 

.00340 

Onloadine  Phase 

16 

22,330 

.00521 

.00507 

.00436 

.00402 

.00280 

17 

14,890 

.00364 

.00354 

.00310 

.00289 

.00221 

18 

7,440 

.00206 

.00202 

.00184 

.00176 

.00161 

19 

0 

.00048 

.00050 

.00058 

.00063 

.00100 

0>) 

Maxlaum 

Applied  Stress  ■ 

35,100  psl 

Loadlna  Phase 

1 

15,190 

.00322 

.00311 

.00257 

.00231 

.00122 

2 

20,050 

.00433 

.00418 

.00349 

.00315 

.00181 

3 

24,920 

.00551 

.00533 

.00450 

.00408 

.00253 

4 

29,780 

.00679 

.00659 

.00652 

.00515 

.00340 

3 

35,100. 

.00836 

.00814 

.00705 

.00652 

.00459 

Unloadlnn  Phase 

12 

26,330 

.00651 

.00634 

.00556 

.00519 

.00388 

13 

17,550 

.00465 

.00455 

.00407 

.00385 

.00318 

14 

8,770 

.00279 

.00276 

.00259 

.00252 

.00247 

15 

0 

.00093 

.00096 

.00110 

.00119 

.00177 

(c) 

Maxliauai 

Applied  Stress  • 

43,200  psl 

Loadlne  Phase 

1 

15,190 

,00322 

.00311 

.00257 

.00231 

.00122 

2 

20,050 

.00433 

.00418 

.00349 

.00315 

.00181 

3 

24,920 

.00551 

.00533 

.00450 

.00408 

.00253 

4 

29,780 

.00679 

.00659 

.00561 

.00515 

.00340 

5 

35,100 

.00836 

.00814 

.00705 

.00652  . 

.00459 

6 

39,150 

.00978 

.00953 

.00836 

.00780 

.00574 

7 

43,200 

.01150 

.01120 

.00999 

.00934 

.00716 

- 

Unloadlnx  Phase 

• 

32,400 

.00890 

.00873 

.00789 

.00748 

.00611 

9 

21,600 

.00661 

.00651 

.00605 

.00584 

.00523 

10 

10,800 

.00432 

.00431 

.00423 

.00420 

.00436 

11 

0 

.00187 

.00193 

.00219 

.00234 

.00307 

COD  at  Center  of  Crack,  in. 
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overly  compressive  effect  of  the  superdislocation  on  the  crack  flanks  that  Is  revealed 
as  the  applied  stress  Is  diminished  which  invalidates  this  computational  approach. 
(Symbols  denote  finite  element  results,  solid  lines  represent  model  solution.) 
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It  was  seen  that,  while  little  discrepancy  existed  in  the  beginning  stages  of  the 
unloading,  this  eventually  became  untenable.  In  particular,  at  zero  load  the 
compressive  effect  of  the  superdislocation  on  the  crack  flanks  was  far. too  great; 
nb,  the  crack-tip  displacement  is  about  right.  Clearly,  to  reduce  this  effect, 
the  superdislocation  strength  must  be  decreased,  or  it  must  be  pushed  further 
away  from  the  crack,  or  both. 

The  set  of  calculations  shown  in  Figure  4 demonstrated  the  Inadequacy 
of  reducing  the  dislocation  strength  while  keeping  its  position  fixed.  For 
these  results  the  strength  was  reduced  by  a factor  of  two,  as  would  be  given  by 
the  procedure  used  previously  for  a zero  load  level.  The  position  was  then 
Increased  to  determine  if  a combination  change  would  suffice.  It  was  found 
that  agreement  with  the  finite  element  results  could  be  obtained  for  redlcal 
strength  changes,  but  only  with  highly  unrealistic  changes  in  the  superdisloca- 
tlon  position. 

The  results  shown  in  Figure  4 indicated  that  focusing  on  the  super- 
dislocation  position  is  the  only  plausible  alternative.  There  is  a physical 
rationale  for  this  which  stems  from  the  point  of  view  that  the  superdislocation 
is  an  idealization  of  a dislocation  pile-up.  Specifically,  its  position  repre- 
sents the  "center  of  gravity”  of  such  a pile-up.  In  such  a pile-up,  the 
dislocations  near  the  crack  will  slip  backwards  and  will  vanish  into  the  crack 
as  the  applied  stress  is  reduced.  In  contrast,  the  dislocations  in  the  pile-up 
away  from  the  crack  will  be  locked  into  position.  Thus,  the  "center  of  gravity" 
of  the  pile-up  will  shift  as  the  load  diminishes.  Equivalently,  the  superdis- 
location representing  the  pile-up  will  shift  outwards.  This,  of  course,  must 
be  accompanied  by  a modest  decrease  in  the  strength. 

While  this  argument  is  appealing,  it  unfortunately  cannot  easily  be 
made  quantitative.  An  alternative  is  to  use  the  finite  element  results  to  learn 
Just  how  the  superdislocation  must  be  shifted.  The  result  of  a trial  and  error 
calculation  demonstrated  that  a monotonlc  Increase  in  the  distance  between  the 
superdislocation  position  and  the  crack  tip,  as  the  load  is  decreased  to  zero 
gives  excellent  agreement  with  the  finite  element  results.  Denoting  the  maximum 
applied  stress  as  and  the  superdislocation  position  at  that  stress  as 

the  positions  I required  for  good  agreement  with  this  set  of  finite  element 

crack-opening  displacements  for  a maximum  load  of  29,780  psl  are  estimated  in 

* 

Table  3.  The  crack-opening  displacement  comparison  resulting  from  movements 

* Note  that  these  results  are  valid  only  for  the  particular  geometry,  load  level,  | 

and  mechanical  properties  used  to  generate  the  finite  element  results  of  Table  2;  | 

NB,  the  fairly  large  aspect  ratio  used  in  these  calculations. 
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ARBITRARILY  VARIED  FOR  HEURISTIC  PURPOSES.  Note  that  position  must  be  increased 
to  an  unreasonable  extent  to  achieve  agreement  with  finite-element  results 
which  invalidates  this  con^jutational  approach.  (Symbols  denote  finite  element 
solutions,  solid  lines  represent  superdislocation  model  solution.) 
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of  the  superdislocations  In  accord  with  Table  3 are  shown  In  Figure  5.  The 
results  shown  In  Figure  5 can  be  readily  seen  to  offer  a significant  Inprove- 
ment  over  those  of  Figure  3. 


TABLE  3.  RELATIVE  SUPERDISLOCATION  POSITION  IN  THE 
UNLOADING  PHASE  OF  A LOAD-UNLOAD  SEQUENCE 
REQUIRED  FOR  CORRECT  PREDICTION  OF  THE  CRACK- 
OPENING DISPLACEMENTS  (See  Figure  5) 


R “ 0 . /a 

min  max 

i/i 

o 

1.00 

1.0 

0.75 

1.5 

0.50 

3.0 

0.25 

10.0 

0 

13.0 

The  approach  to  be  taken  In  the  following  section  is  based  on  the 
point  of  view  that  the  residual  plasticity  created  during  an  extended  period 
of  steady-state  crack  growth  under  constant  AK  load  cycles  can  be  represented 
by  a single  superdlslocatlon  pair.  The  strength  of  the  "residual"  superdls- 
locatlons  Is  governed  by  Equations  (8)  and  (9) . Its  position  Initially  arises 
from  these  same  equations  but  Is  modified  to  match  the  finite-element  solutions 
for  the  crack  faces  during  unloading.  An  aid  In  obtaining  these  results  Is 
the  linearized  relation  for  the  crack-opening  displacement  at  the  center  of  the 
crack.  This  Is  obtained  by  using  the  relation 

d(l.O)  - 2 - I (ff  cos  I . (17) 

where  d ~ d(l,x)  Is  the  one-superdislocatlon  version  of  the  function  of  a com- 
plex variable  Included  In  Equation  (7).  Substituting  Equation  (17)  Into  (7)  then 
gives 
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Note  that  v represents  here  the  total  distance  between  the  crack  faces  and 
Is  twice  the  distance  from  the  crack  faces  to  the  symmetry  plane. 
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4.  FATIGUE  CRACK-GROWTH  CALCULATIONS 


In  previous  work,  fatigue  crack  growth  computations  were  carried  out 
using  one  superdlslocatlon  pair  to  represent  the  crack-tip  plasticity  In  each 
and  every  load  cycle.  This  can  be  referred  to  as  the  "cycle-by-cycle"  model. 

In  contrast,  the  following  calculations  will  be  performed  using  a much  simpler 
two  superdlslocatlon  model.  This  Is  motivated  by  the  necessity  of  working  with 
a less  cumbersome  coiq;>utatlonal  scheme  In  the  formative  stages  of  the  work  and 
the  fact  that.  In  the  cycle-by-cycle  model,  one  of  the  superdlslocatlons  left 
behind  by  the  growing  crack  Invariably  dominates.  This  suggests  that  a suitable 
model  for  steady-state  crack  growth  under  constant  AK  load  cycles  would  be  one 
composed  of  one  "residual"  superdlslocatlon  pair  and  a nascent  superdlslocatlon 
pair.  This  model  Is  Illustrated  In  Figure  6. 

Formulation  of  the  Simplified  Model 

The  relations  given  earlier  apply  to  the  single  superdlslocatlon  pair 
representation  of  residual  crack-tip  plasticity.  For  the  steady-state  crack- 
growth  situation  (which  will  occupy  the  balance  of  this  section).  It  Is  the 
Interaction  of  a nascent  superdlslocatlon  pair  with  the  residual  plasticity 
that  Is  of  chief  concern.  The  relations  for  the  Interaction  between  two  super- 
dlslocatlons are  also  necessary,  therefore.  Linearized  versions  of  the  dislo- 
cation Interaction  terms  can  be  deduced  by  assuming  that  both  dislocations  are 
on  the  same  slip  plane  and  that  they  occupy  positions  such  that  a >> 

Hence,  the  Interaction  terms  are  found  to  be 


sln^  6 cos^  -I-  (19) 


and 


i * i, 

*21  ^ *^21  “ ^ ® 1 ^ — \ ~ ^ 

1^2  *■! 


(20) 


from  which  It  can  be  seen  that 

*21  ‘'21  “ ■ ^*12  *'12^ 
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FIGURE  6.  SIMPLIFIED  TWO  SUPERDISLOCATICW 
MODEL  FDR  STEADY-STATE  FATIGUE 
CRACK  GROWTH 
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For  convenience  In  the  following,  the  specific  values  for  all  the  linearized 
functions  for  0 * cos  ^ (1/3)  can  be  written  out.  These  are 


fi  = 6.532  (o/fi)^ 

hji  = 0.2722  + 

8i3  '’iz 


fs  - 6.532  (a/fj)^ 

8s2  + ^22  = “ 1/^2 

0.314  hj  - 0.2722  (c/Zg)^ 

- 1 128  _l_ 

■ ■ ^ (£,^2)^  ■ 81  X, 


+ 0.314 


821  1^21 


16  (£i£a)^  128  £2 

9 £1=  81 


A comparison  showing  the  accuracy  of  the  latter  two  functions  Is  shown  In 
Reference  [7]. 

Just  as  In  the  previous  work,  the  equations  for  the  model  governing 
the  conditions  at  maximum  load  are  for  the  force  equilibrium  of  each  super- 
dislocation  and  for  the  cancellation  of  the  singularity.  The  latter  Implies 
that  crack  growth  has  occurred  on  the  Increasing  load  part  of  the  load  cycle 
and  that,  when  the  growth  ceases,  the  crack-drlvlng  force  Is  just  zero.  For 
this  reason.  It  Is  Inappropriate  to  Impose  the  singularity  cancelling  equation 
at  the  minimum  load  as  was  done  previously.  Hence,  this  Is  no  longer  done. 
Another  major  change  In  the  model  formulation  Is  to  allow  a repositioning  of 
residual  superdislocations  at  the  minimum  load  as  discussed  above.  This  Is  based 
on  the  Idea  that  a superdislocation  represents  the  "center  of  gravity"  of  a 
dislocation  plle-up.  In  such  a plle-up,  the  dislocations  near  the  crack  will 
slip  backwards,  vanishing  Into  the  crack  as  the  applied  stress  Is  reduced.  The 
remaining  dislocations  cannot  overcome  the  friction  stress  and  are  locked  In 
position.  Thus,  using  the  center  of  gravity  point  of  view,  the  superdislocation 
representing  the  residual  plasticity  at  the  minimum  load  must  be  further  from 
the  crack  plane  and  will  have  a somewhat  reduced  strength. 

To  obtain  a solution  for  the  steady-state  situation,  consider  that 
the  residual  superdislocation  Is  located  at  a position  L,  the  nascent  superdis- 
location at  a position  £.  A solution  exists  if,  and  only  if,  the  nascent 
superdlslocatlon  Is  In  force  equilibrium  and  the  net  stress-intensity  factor 
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Is  zero.  For  the  latter  condition  to  hold  with  a nonnegative  value  of  the 
strength  of  the  nascent  superdislocation,  there  will  be  a threshold  value  of 
the  applied  stress.  From  Equation  (5),  for  the  residual  superdislocation 
strength  held  constant,  thla  Is 


a 


th 


(21) 


where  L/L  Is  the  ratio  of  the  adjusted  position  to  the  Initial  position  of 
o 

the  residual  superdislocation.  Now  writing 


<“>th  - - ".rin’ <•“)'’ 

It  can  be  seen  that 


(22) 


(23) 


where,  as  usual,  R * K 

mn  iDAx 

A "threshold"  will  exist  in  the  model  If  the  right-hand  side  of 
Equation  (23)  Is  positive.  This  will  be  the  case  If 


Using  the  values  of  L/L^  given  in  Table  3,  Figure  7 shows  that  the  Inequality 
Is  always  satisfied  for  all  R < 1.  Thus,  the  threshold  effect  does  exist  In 
the  model. 
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O Position  adjusted  in  accord! 
with  the  finite  element 
COD  results  - see  Table  3 
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FIGURE  7.  COMPARISON  OF  THE  RESIDUAL  SUPERDISLOCATION 
POSITION  WITH  VALUES  FOR  ZERO  THRESHOLD 
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It  might  be  noted  that  this  arises  strictly  from  consideration  of  the 
residual  compressive  stresses  exerted  on  the  crack  tip  by  the  plasticity  generated 
In  previous  load  cycles  and  Is  Independent  of  a crack-closure  condition.  It  can 
be  further  seen  that  there  is  nothing  In  the  nature  of  a material  property  as- 
sociated with  the  threshold  as  some  researchers  have  suggested. 

The  final  step  in  the  formulation  of  the  steady-state  solution  Is  to 
determine  b.  Then,  by  associating  the  crack-advance  Increment  with  the  length 
of  new  crack  surface  projected  along  the  original  crack  plane 

Aa  ” 2k b cos  6 , (24) 

where  k Is  taken  to  be  an  environment-material  constant.  The  parameter  k Is 
Intended  to  account  for  the  experimentally  observed  fact  that  "rewelding"  of 
the  crack  surfaces  can  occur  and,  further,  that  the  extent  to  which  this  occurs 
depends  on  the  environment.  Thus,  0 < k _<  1.  In  benign  evlronments,  k would  be 
expected  to  be  quite  small.  In  agresslve  environments.  It  would  approach  unity. 

Steady-State  Crack-Growth  Rate 

Using  the  simplified  two  dislocation  model  of  Figure  6,  a prediction 
of  steady-state  fatigue  crack  growth  can  be  made.  For  deflnltness,  the  follow- 
ing parameters  will  be  used: 

B - strength  of  the  residual  superdislocation 
b - strength  of  the  nascent  superdislocation 

L • position  of  the  residual  superdislocation  along  the  slip  plane 
I - position  of  the  superdislocation  along  the  slip  plane 
c * distance  from  the  crack  tip  to  the  residual  superdislocation 
slip  plane. 

It  will  be  considered  that  B,  L,  and  c are  fixed  by  the  prior  load  history;  e.g.. 
Independently  by  finite  element  or  other  computations.  The  variables  to  be  de- 
termined are  b and  1.  This  Is  done  by  satisfying  the  equations  for  the  equilibrium 
of  the  nascent  superdislocation  and  for  the  cancellation  of  the  singularity  at 
the  maximum  load  In  the  cycle.  Thus,  from  Equations  (8)  and  (9) 

\ " °max  ka,^)|+  glg(L,t)  + k (L,Mi  (25) 
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and 


aa^  {b£(a)  + Bf(L)}  . (26) 

UdX  wTT 

The  parameter  b can  be  determined  from  Equation  (26)  as 

Then,  by  using  Equatlon(27)  to  eliminate  b from  Equation  (25),  I can  be  determined. 
The  result.  In  terms  of  the  linearized  versions  of  the  functions  f,  h,  g,  and  k 
as  given  by  Equations  (10-13) , Is 


1 

( 1 / 

, r 

_1 

1 

( ^ ^ 

2/  max  1^ 

,5  S xTi  w COS  2 ^ 

sin  0 cos 

1 EB  r 
8it  ^ L 


1-3  sln^  6 cos^ 


. - . . 2 (1  - cos  4 0 ) EB  i 

T . - a sin  e cos  0 + 5 -j-  > 

I 1 max  Bn  L I 


-2 


(28) 


The  equivalent  linearized  version  of  Equation  (27)  is 


E ) 6 


/a\**  '^max 

1 

Uj  ^ „ 0 

8ir 

A 

' ' sin  0 cos 

(29) 


Combining  Equations  (28)  and  (29)  then  gives 

” - A p - -TT — Tol!  Vx  -{iif  ® I ^ 

L ain  0 cos  -rjl  L 


. . . . 1 - cos  4 0 EB 

• ^ ^1  - ®max  6 ® + 4^1^ T 


-1 


(30) 
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Setting  K « a (ir*)^  and  t.  » Y sin  6 cos  6 then  gives 
® max  max  i 


9EY  sin  6 cos  6 


[1  /l  - cos  4 8 \ EB 1 ^ 

^ ” Y Air  \sln  0 cos  0/  YLJ  ’ 


(31) 


where 


“ (ife)  S ® ‘'°®  I 


th 


(32) 


Further  progress  can  be  made  by  relating  B and  L to  the  values  for  fixed  load 

conditions.  Call  these  B and  L , respectively.  It  is  readily  found  that 

o o 


B 


max 


o 9EY  sin  0 cos  0 


(33) 


and 


K Q 1 

max  „ 0 ] 

L “ 2 2 2 — ^ ® 7 

° 72irY  sin  0 cos  0 sin  0 


2 

a 1 

1 

. ■ M 

-2 


(34) 


Key  parameters  are  seen  to  be  the  ratios  B/B^  and  h/L^.  These  can  be  conveniently 
expressed  in  terms  of  dimensionless  parameters  B and  defined  as 


B 


■(t) 


and 


which  fix  the  residual  plasticity  In  this  model. 
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Before  specializing  the  results  for  0 » cos  ^(1/3),  It  might  be 
emphasized  that.  In  order  for  the  model  to  correctly  predict  general  yielding 
when  the  applied  stress  Is  equal  to  the  tensile  yield  stress  Y,  the  friction 
stress  was  taken  equal  to  Y sin  6 cos  6 • 0.314  Y.  Further  note  that  plane 
stress  conditions  are  used  and  that  the  equations  have  been  linearized;  l.e., 
by  taking  1 <<  L <<  a,  and  neglecting  the  distance  between  the  two  slip  planes. 
Under  these  conditions,  substituting  Equation  (31)  Into  Equation  (24)  then  gives 


da  _ 7 

dN  22^ 


EY| 


_ K 

j.  ®max  . 4^2’  EB 

1 max 

th 

t P " Y 9tt  YL 

-1 


(37) 


The  Introduction  of  expressions  for  B and  L can  be  facilitated  by  expressing 
them  in  terms  of  the  values  obtained  for  a single  load  application  through  the 
use  of  Equations  (33-36).  In  particular 


B - 


and 


L 


147  ^ max 

1024Tr  ^2 


(38) 


(39) 


where  B and  A,  as  above,  denote  the  ratios  of  B and  L for  a growing  crack  (l.e., 
one  with  a wake  of  plasticity)  to  their  values  for  a single  load  application. 

Using  Equations  (38)  and  (39)  to  eliminate  b and  L In  Equation  (37) 
then  gives  a result  that  can  be  written  as 


da 

dN 


K 

max 


» 


(40) 
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f 


\ 

\ 


f 

I 


where 


Y s ^ Y 

^ max 


and 


ZjLji  + 

32>15‘1 


-1 


(41) 


(42) 


Because  K - K “ (K  - K ) - (K  “ K j ) = K - (AK)_,  for  small  applied 

stresses.  Equation  (40)  can  be  seen  to  be  qualitatively  consistent  with  Equation 

(1).  Quantitative  comparisons  will  depend  on  the  exact  form  of  the  functions 

6 and  X.  While  these  are  expected  to  be  principally  dependent  on  R and  on  Y/E, 

general  relations  have  not  yet  been  determined.  However,  It  Is  Interesting  to 

note  that,  if  B/X  is  on  the  order  of  unity  and  k <•  1,  then,  from  Equation  (42), 

A will  be  approximately  equal  to  .05,  a value  In  quite  good  agreement  with 

McEvlly’s  findings  for  crack  growth  In  air. 

There  are  some  Important  differences  between  Equation  (1)  and  (40) 

also.  First,  In  this  work,  K . emerges  as  a variable  which  reflects  a load 

th 

history  effect.  In  contrast,  the  conventional  view  of  the  threshold  is  that  it 
Is  a material  property.  Determinations  of  this  quantity  are  made  In  such  a way 
as  to  minimize  the  load  history  effect  by  gradually  diminishing  the  stress 
amplitude  to  determine  the  limit  of  zero  crack  growth.  Clearly,  the  value  given 
by  Equation  (41)  Is  appropriate  to  a particular  load  history,  not  to  this 
limiting  case. 

A second  difference  Is  that,  while  a term  does  not  appear  In  the 
model,  an  explicit  dependence  on  does  emerge  from  this  approach.*  It  Is 

possible  that  another  crack  growth  mode  cond>lnes  with  the  sliding  off  mode 
being  considered  here  at  high  values.  But,  It  Is  difficult  to  see  how 

this  could  be  connected  with  a value  In  a logical  self-consistent  manner. 

* As  shown  in  Reference  [11],  Equation  (1)  can  be  generalized  to  take  account 
of  accelerated  crack  growth  as  net  section  yielding  Is  approached. 
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The  Crack  Closure  Condition 

An  essential  difference  between  the  previous  work  ard  the  present  is 
with  regard  to  the  crack-closure  prediction  capability  of  the  model.  It  was 
stated  in  Reference  7 that  the  imposition  of  the  singularity  cancelling  equation 
forces  the  crack  faces  to  close  smoothly  at  the  tip  and,  consequently,  that  this 
is  a necessary  condition  for  crack  closure  to  occur.  But,  as  the  work  given 
in  Appendix  A shows,  the  interpretation  of  the  crack-opening  displacement 
relations  giving  rise  to  this  conclusion  was  faulty. 

It  is  now  clear  that  a step  change  (of  height  equal  to  the  superdis- 

it 

location  strength)  occurs  in  the  crack-opening  displacement.  Thus,  while  can- 
celling the  crack- tip  singularity  still  forces  the  crack  faces  to  be  parallel, 
they  will  not  touch  each  other  except  in  the  limiting  case  where  the  superdis- 
location strength  is  zero.  It  can  be  shown,  in  fact,  that  no  impingement  of 
the  surfaces  is  possible  at  any  point  along  the  crack  flanks  for  a positive 
superdislocation  strength.  Proof  of  this  statement  is  given  in  Appendix  A of 
this  report.  However,  while  it  may  appear  that  crack  closure  is  impossible  in 
this  model,  that  this  is  not  the  case  can  be  made  clear,  as  follows. 

For  a superdislocation  on  a slip  plane  that  intersects  the  crack  faces 
at  the  crack  tip,  no  crack  closure  is  possible.  However,  the  possibility  of 
crack  closure  does  exist  if  the  slip  plane  intersects  the  crack  faces  at  any 
other  point.  This  is  Illustrated  in  Figure  8.  This  can  be  demonstrated 
quantitatively  by  use  of  the  linearized  version  of  the  crack-opening  displacement 
at  the  point  of  Intersection.  The  key  is  the  relation 

d(L,a-c)  - - ) cos  I (43) 

where  c denotes  the  distance  between  the  crack  tip  and  the  point  of  intersection 
of  the  superdislocation  slip  plane  and  the  crack  face;  l.e.,  the  point  x - a - c. 

Note  that  Equation  (43)  is  based  on  the  assumption  that  c < L as  well  as  that 
L <<  a. 

* The  relations  exhibited  in  References  [7-9]  are  all  correct  provided  the  Heaviside 
step  function  H(x)  replaces  the  Dirac  delta  function  6(x)lnthe  crack-opening 
displacement  equations. 
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figure  8.  SUPERDISLOCATION  MODEL  OF  RESIDUAL  CRACK- TIP 

PLASTICITY  SHOWING  THE  POSSIBILITY  OF  CRACK  CLOSURE 
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) 


Substituting  Equation  (43)  Into  Equation  (7)  will  give  an  expression 
for  the  crack-opening  displacement  at  the  potential  point  of  contact  between 
the  two  crack  surfaces.  For  plane  stress  conditions  (which  are  appropriate 
because  closure  effects  are  usually  observed  on  the  surface) , this  Is 

Y V » (2ac)’  - sin  e cos  j , (44) 


where  the  elastic  contribution  has  also  been  llnearlezed  for  consistency.  Using 
Equation  (44),  for  a given  B and  L,  the  applied  stress  at  which  contact  of  the 
crack  faces  will  occur  can  be  readily  determined.  This  Is 


a 

c 


_3_ 

2ii 


EB 

(2a  L)** 


sin  0 cos  Y 


(45) 


Note  that,  although  Equation  (45)  Is  apparently  Independent  of  the  yield  stress 
and  of  the  distance  c,  this  Is  misleading  because  B may  possibly  be  a function 
of  these  parameters . 

Defining  • a^(7ra)  and  using  Equation  (45)  to  express  a^.  It  Is 
found  that  ' 


“ (ii)  ^ ® 2 

L 

It  can  be  seen  that  this  result  Is  Just  the  same  as  the  threshold  value  given 
by  Equation  (32).  Thus,  It  would  seem  that  whether  one  takes  the  point  of  view 
that  residual  stress  or  crack-closure  effects  govern  fatigue  crack-growth  rates 
may  be  Irrelevant — both  effects  may  occur  simultaneously.  It  Is  also  worth 
repeating  that  both  quantities  are  simply  reflections  of  the  residual  plasticity 
In  the  material  and  are  not  associated  with  any  material  properties  other  than 
the  usual  mechanical  parameters. 

To  reiterate,  while  the  crack-closure  condition  arising  previously  In 
the  model  Is  now  Inoperative,  the  phenomenon  Is  contained  In  the  present  model 
In  another  way.  Again,  this  occurs  In  the  model  In  a natural  way  and  Is  not 
forced  upon  it  as  a basic  postulate. 
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Comparison  with  Experimental  Results 

The  strength  and  position  of  the  residual  superdislocation  obviously 
exert  strong  Influences  on  the  crack-growth  rate.  It  Is  believed  that,  by 
extending  the  procedure  described  In  Section  3 of  this  report,  the  values  of 
these  parameters  can  be  determined  by  comparison  with  elastic-plastic  finite 
element  solutions.  Specifically,  a series  of  load,  crack  growth,  unload  compu- 
tations can  be  performed  In  which  growth  progresses  (by  releasing  nodes)  with 
load  shedding  to  achieve  a constant  K value.  Values  of  $ and  A can  then  be 
determined  by  matching  the  crack-face  displacements  during  the  unloading.  This 
has  not  yet  been  accomplished,  however. 

As  an  e^edlent,  to  test  the  results  of  the  model  against  experimental 
results  on  constant  aiiq>lltude  fatigue  crack  growth,  an  alternative  approach 
will  be  used  to  set  values  of  3 and  A.  Guided  by  the  results  described  In 
Section  3,  for  simplicity.  It  can  be  assumed  that  3 ~ 1.  The  well  known 
experimental  result  of  Elber  [12]  can  then  be  used  to  evaluate  A.  Specifically, 
Elber  found  for  2024-T3  aluminum  alloy  that 


max 


- K 


_op_ 


max 


- K 


0.5  + 0.4R 


‘min 


or 


op 


(0.5  + O.IR  + 0.4R  ) K 


max 


(47) 


It  Is  also  possible  to  determine  K from  the  model. 

op 


As  shown  In  the 

foregoing,  this  Is  done  by  determining  the  load  level  at  which  the  COD  at  the 
point  of  Intersection  of  the  residual  slip  plane  with  the  crack  Is  just  equal 
to  zero;  cf.  Equation  (46).  Using  Equations  (38)  and  (39)  then  gives 


K 


op  max 


(48) 


Combining  Equations  (47)  and  (48),  setting  3 ~ 1,  gives 


A - (0.5  + O.IR  + 0.4R^)'‘^ 


(49) 
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Finally,  substituting  Equation  (49)  Into  Equation  (40),  making  use  of  (41)  and 
(42),  then  gives  the  prediction 


dN 


0.155 

EY 


r 0.5-.0.,.lR-0..4  1 

Li+2.03  (0.5  + 0.1R+0.4R^)^ J 


K‘ 


max 

TlTo  TyT 

max 


(50) 


for  the  steady-state  crack  growth  rate  resulting  from  constant  amplitude  loading. 

A particularly  well  documented  set  of  experimental  data  that  can  be  used 

to  test  Equation  (50)  Is  that  of  von  Euw,  et  al  [13]  on  2024-T3  aluminum.  Their 

data  were  obtained  from  a 3-lnch  wide  SEN  specimen  at  two  different  crack  length 

to  specimen  width  ratios,  a/w.  As  they  pointed  out,  the  conventional  view  Is  that, 

for  the  same  values  of  K and  K . , the  same  crack  growth  should  be  obtained 

1DAX  mn 

for  all  a/w  values.  However,  this  was  not  the  case.  To  determine  If  Equation 

(50)  could  account  for  their  observed  differences,  the  (1  - ^ terms  was 

Interpreted  as  a net  section  value  with  a calculated  from  the  given  K 

max  max 

values  via 


, -K  (i.s)£(2)r‘  , 

max  maxi  Vw/  ' w'  'w') 

where  f(a/w)  Is  the  finite  width  correction  for  the  SEN  specimen  used  In 
Reference  [13].  In  particular,  f(0.25)  " 2.308  and  f(0.45)  “ 4.989  which  gives 

o ■ K /78.7  and  a ■ K /167.4  for  a/w  - 0.25  and  0.45,  respectively. 

max  max  max  max  ^ 

Other  values  used  In  the  calculation  are  E * 10  ksl  and  Y • 52.5  ksi,  also  as 
given  In  Reference  [13]. 

Comparisons  between  the  prediction  of  Equation  (50)  and  the  results  of 
von  Euw,  et  al  for  two  different  kinds  of  tests  are  shown  In  Figures  9 and  10. 

In  Figure  9 R Is  a constant  while  Figure  10  shows  a case  where  AK  Is  constant. 

In  both  comparisons  It  can  be  seen  that  the  model  does  an  excellent  Job  of 
predicting  the  growth  rates.  Also,  the  previously  unexplained  difference  be- 
tween the  results  for  the  two  different  a/w  ratios  Is  now  accounted  for. 

Note  that  the  data  used  in  Figures  9 and  10  are  Just  the  same  as  those 
used  to  test  the  earlier  version  of  the  model — see  Figures  5.8  and  5.9  (page  56 
snd  57)  of  Reference  [7].  In  fact,  the  discrepancies  obtained  were  a primary 
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motivation  for  the  further  development  that  has  led  to  the  present  model. 

Clearly,  the  model  is  now  in  complete  accord  with  the  well-established  experi- 
mental trends  exemplified  by  the  data  of  von  Euw,  et  al. 

It  should  perhaps  be  emphasized  that  the  development  given  in  this 
section  applies  only  to  constant  amplitude  loading.  While  the  model  can  certainly 
be  applied  to  variable  amplitude  loading — and  this  is  the  ultimate  goal  of  the 
work — such  computations  will  require  that  additional  degrees  of  freedom  be 
^ allowed;  l.e.,  more  than  just  one  residual  superdislocation  pair.  It  is  now 

clear  how  this  can  be  done.  As  for  example,  as  described  in  pages  58-62  of 
Reference  [7],  heuristic  calculations  Indicate  that  meaningful  results  will  be 
obtained.  However,  more  detailed  calculations  such  as  these  are  beyond  the  scope 
of  this  work. 


! 
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5.  Discussion.  CONCLUSIONS.  AND  REC(»<MEtn)ATIOMS 

The  work  described  herein  was  directed  at  the  satisfactory  completion 
of  the  following  major  goals: 

(1)  Prediction  of  constant  crack-growth  rates  for  load  levels 
that  are  uniform  from  cycle  to  cycle  that  are  In  satisfactory 
agreement  with  experimental  results. 

(2)  Predictions  that  properly  reflect  the  effect  of  mechanical 
properties  on  fatigue  crack-growth  rates. 

(3)  Computational  efficiency  such  that  a large  number  of  load 
cycles  can  be  performed  In  a reasonable  amount  of  time. 

In  addition,  as  has  been  the  Intent  from  the  beginning  of  the  work,  the  mathe- 
matical model  Incorporates  only  material  properties  for  which  data  can  readily 
be  found  for  the  materials  of  Interest  In  engineering  applications.  While  the 
model  Is  not  yet  complete  In  all  respects,  the  results  provided  In  this  report 
demonstrate  that  these  goals  have  essentially  been  met. 

One  criticism  of  the  results  presented  In  Reference  [7]  was  that  the 
crack-growth  rate  always  turned  out  to  be  Inversely  proportional  to  the  yield 
stress.  While  there  Is  no  decisive  body  of  experimental  evidence  which  reveals 
the  proper  dependence  on  Y,  there  are  some  data  which  Indicate  that  a Y de- 
pendence Is  not  universal.  The  present  results  partly  alleviate  this  deficiency 
in  that  the  yield  stress  dependence  now  differs  from  the  Y ^ types  of  dependence: 
see,  for  example.  Equation  (40).  Further  Improvements  are  also  possible,  as 
follows . 

The  current  model  represents  crack-tip  plasticity  only  In  terms  of  the 
yield  stress  of  the  material.  This  means  that  it  is  confined  to  elastlc-perfectly 
plastic  behavior.  More  realistic  material  behavior  might  be  Incorporated  Into 
the  model  by  following  a suggestion  made  by  Rice  [14].  This  Is,  for  a power  law 
hardening  material,  l.e.,  stress  proportional  to  strain  to  the  N^  power  In  the 
plastic  regime,  the  tensile  flow  strength  corresponding  to  the  Mlses  equivalent 
shear  strain  N/(l+N)  can  be  used.  This  Is 
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where  Y Is  the  yield  stress,  E Is  the  elastic  modulus,  and  v Is  Poisson's  ratio. 
A more  convenient  form  Is  obtained  by  using  the  shear  modulus  G.  This  Is 


o 

flow 


gN 


(N 

"^l.  Consequently,  the  latter  form  clearly  shows  that 
for  a nonstrain  hardening  material,  “ Y.  The  Introduction  of  such  a 

relation  Into  the  model  would  require  four  basic  mechanical  properties  to  specify 
material  behavior:  G,  v,  Y,  N. 

The  above  provides  one  means  by  which  more  realistic  elastic-plastic 
behavior  could  be  Incorporated  Into  the  model.  Another  Is  to  use  cyclic  stress- 
strain  properties  rather  than  the  monotonlc  values.  In  terms  of  the  single 
two-superdlslocatlon  model  developed  In  the  above,  the  same  value  of  need 
not  be  used  for  each  superdlslocatlon.  For  exan^le,  the  residual  superdlsloca- 
tlon  can  be  associated  with  the  uniaxial  tensile  yield  stress,  while  the  nascent 
superdlslocatlon  can  be  associated  with  the  cyclic  value.  This  may  produce  a 
delayed  minimum  In  the  crack-growth  retardation  following  an  overload  as  the 
location  of  the  nascent  dislocation  changes  from  virgin  material  to  cyclically 
deformed  material.  There  are  of  course  other  explanations  for  a delayed  minimum; 
e.g. , difference  between  crack  closure  at  the  surface  and  In  the  Interior.  The 
point  here  Is  that  such  freedom  exists  and  can  be  exploited  In  subsequent  stages 
of  development  of  the  model. 

A number  of  other  research  efforts  currently  In  progress  are  also 
aimed  at  developing  a general  fatigue  crack-growth  model.  These  all  utilize 
either  a sl]iq>llfled  representation  of  the  crack-tip  plasticity  or  Incorporate 
a mechanism  for  the  crack-closure  condition  or  both.  Evans  [15]  uses  a simplified 
version  of  the  superdlslocatlon  concept  In  his  model  for  stress  corrosion  crack- 
ing and  fatigue.  Paris  and  Herman  [16]  have  developed  an  approximate  disloca- 
tion baaed  model  that  predicts  crack-closure  and,  hence,  an  effective  AK  that 
appears  to  be  In  good  accord  with  experimental  results.  Newman  [17]  has  devised 
a ligament  model  based  on  finite  element  calculations  [18]  which  leaves  plastically 
deformed  material  In  the  wake  of  the  advancing  crack.  Most  recently,  Budlansky 
and  Hutchinson  [19]  have  extended  a Dugdale  model  crack-tip  plasticity  representa- 
tion to  calculate  crack-clostire  effects  In  steady-state  crack  growth. 
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The  essential  difference  between  the  work  described  In  this  report 
and  the  work  of  others  Is  that  the  latter  are  not  able  to  Incorporate  a crack- 
growth  criterion  directly  Into  their  models.  The  model  described  here  appears 
to  be  uniquely  capable  of  providing  this  key  feature.  As  a consequence,  the 
reliance  on  a postulated  controlling  mechanism  (e.g.,  crack-closure,  residual 
plasticity)  otherwise  necessary  Is  avoided.  These  mechanisms  are  automatically 
Included  In  this  model  and  make  their  contribution  In  a natural  way.  It  might 
further  be  noted  that  experimental  observations  of  fatigue  crack  tips  also 
provide  support  for  the  formulation  of  the  model.  The  mechanisms  observed  by 
Lankford  and  Kusenberger  [20]  and  by  Neumann  [21]  can  be  cited  In  this  regard. 

In  conclusion.  In  this  report,  two  areas  In  which  the  formulation  of 
the  model  were  Inadequate  have  been  Identified  and  Improved  upon.  These  Im- 
provements have  led  to  a "steady-state"  fatigue  crack-growth  model  that  Is  In 
good  agreement  with  experimental  data.  Because  these  predictions  can  be  made 
with  a simple  two-superdislocatlon  model.  It  appears  that  the  desired  computational 
convenience  Is  also  at  hand.  However,  further  work  Is  needed  to  qualify  the 
model  for  routine  engineering  calculations.  IVo  specific  recommendations  can 
be  made.  The  first  Is  required  In  order  to  cope  with  variable  amplitude  load- 
ing, the  second  to  cope  with  environmental  effects.  These  are  as  follows. 

Variable  aiiq>litude  loading  will  produce  residual  plasticity  behind 
the  tip  of  a fatigue  crack  that  will  depend  upon  the  details  of  the  load  history. 
The  model  cannot  be  expected  to  reflect  this  with  only  a single  residual  super- 
dislocation  pair  as  was  possible  for  constant  amplitude  loading.  But,  the 
manner  In  which  the  simple  model  can  most  effectively  be  expanded  to  treat 
nonuniform  loads  without  requiring  the  complexity  of  a full  cycle-by-cycle 
representation  of  crack-tip  plasticity  is  not  yet  established.  Further  work  Is 
needed  to  develop  a set  of  rules  for  coiid>lnlng  and  repositioning  superdlsloca- 
tlons  to  represent  the  residual  plasticity  In  hlgh-to-low  and  low-to-hlgh 
sequences.  This  offers  the  possibility  of  treating  multiple  and  block  over- 
loads, negative  excursions,  and  other  complications  Involved  In  service  con- 
ditions. The  authors  believe  that  the  model  described  In  this  report  is  uniquely 
capable  of  providing  this  generality. 

A key  feature  of  fatigue  Is  the  quite  Important  role  played  by  the 
environment.  The  model  developed  here  could  account  for  this  by  taking  the 
Internal  friction  stress  to  depend  on  the  environment.  Alternatively,  a more 
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usefiil  approach  may  be  simply  to  take  some  fraction  of  the  crack-tip  crack-opening 
displacement  as  the  crack-advance  Increment.  The  physical  basis  for  the  latter 
Is  simply  that  the  CTOD  represents  new  surface  area  which  can  be  attacked  by 
the  environment  to  which  It  Is  exposed.  Thus,  "reweldlng"  of  the  surfaces  as 
observed  In  the  experimental  work  of  Neumann,  et  al  [21,22],  can  be  di.rectly 
taken  Into  account.  Benign  environments  producing  little  deterent  to  reweldlng 
will  have  low  growth  ratios  while  detrimental  environments  will  produce  large 
growth  ratios.  Further  work  Is  needed  In  which  the  predictions  of  the  model  are 
tested  by  detailed  comparisons  with  e:q>erlmental  results  to  Implement  and  verify 
this  approach. 
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ADAPTATIONS  OF  FOBMULAS  FOR  THE  STRIP-YIELD 
SUPERDISLOCATION  MODEL  FOR  CRACK-TIP  PLASTICITY 


Introduction 


The  fatigue  crack  propagation  model  described  In  this  report  depends 
on  calculations  from  formulas  based  on  the  strip-yield  superdlslocatlon  model 
for  crack-tip  plasticity  given  In  Reference  [A-1].  The  formulas  are  derived 
as  complex  functions  of  a coaq>lex  variable.  Hence,  computations  from  them  tend 
to  be  slow  and  may  possibly  Involve  obscure  ambiguities.  The  formulas  are  sub- 
jected here  to  rearrangement  Into  real  form  In  the  hope  of  finding  ways  to  verify 
the  results  to  shorten  confutations  based  on  them,  and  to  facilitate  the  develop- 
ment of  simple  approximate  relations. 

From  one  point  of  view,  the  formulas  of  Reference  [A-1]  are  based  on 
the  presumption  of  superdlslocatlons  of  seml-lnflnlte  length  Interacting  with  a 
crack.  To  obtain  correct  results,  the  "tails"  of  the  superdlslocatlons — l.e., 
the  position  of  the  slip  plane  on  the  opposite  side  of  the  crack  plane  from  the 
superdlslocatlon — must  therefore  be  abolished.  There  Is  another  point  of  view, 
however.  This  Is  that  the  previously  derived  relations  can  be  equally  well 
derived  In  terms  of  slip  planes  terminating  In  the  crack  plane.  This  gives  the 
same  stresses,  but  alters  the  displacements.  To  verify  this,  the  following 
proceeds  from  the  first  point  of  view. 


Background 

Consider  that  the  crack  Is  stress  free  while  a load  normal  to  the  crack 
Is  applied  remotely.  All  the  superdlslocatlons  are  Inclined  to  the  crack  by  the 
same  angle  or  Its  supplement,  depending  on  the  quadrant  In  which  the  branch  point 
of  the  superdlslocatlon  lies.  The  entire  configuration  Is  symmetric  about  the 
plane  and  center  of  the  crack.  The  strength  of  each  superdlslocatlon  Is  presumed 
to  be  uniform  throughout  the  entire  length  of  that  superdlslocatlon.  The  nota- 
tion used  for  a single  set  of  superdlslocatlons  Is  shown  In  Figure  A-1. 
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The  equations  of  Interest  furnish  terms  for  a possible  singularity 
cancelling  equation  (for  stress  at  the  crack  tip  considering  effects  of  all  sets 
of  superdislocations)  and  for  determination  of  the  shear  stress  t active  along 
each  superdls location  at  Its  branch  point  (considering  effects  of  the  load  and 
all  sets  of  superdislocations).  In  particular,  with  M sets  of  superdislocations, 
terms  fj,  k^,  k^^  are  to  be  supplied  for  the  conditions 
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The  factor  (l-y)  Is  Introduced  In  case  only  partial  satisfaction  of  the  singularity 
cancelling  equation  Is  desired.  In  addition,  terms  d^  are  to  be  supplied  to 
allow  calculation  of  the  crack-opening  displacement 
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Before  proceeding  further.  It  may  be  observed  that  the  superdislocations 
presumed  In  the  derivation  of  the  formulas  shown  In  Reference  [A-l[  had  seml-lnflnlte 
length  which  is  not  really  desired.  To  make  the  superdislocations  terminate  at 
the  crack  plane,  one  needs  only  to  subtract  effects  contributed  by  the  tails 
beyond  the  crack  plane  (that  Is  the  dotted  parts  of  the  lines).  These  effects 
can  be  found  from  the  existing  formulas  by  letting  1 0 for  Influential  (as  opposed 

if 

to  Influenced)  dislocations.  A second  observation  Is  that  In  their  only  relevant 
context,  the  functions  k^^  and  gj^  are  added,  so  that  they  lack  Individual  signifi- 
cance. Yet  another  observation  Is  that  effects  of  the  superdislocations  at  the 
left  end  of  the  crack  apparently  have  been  omitted  from  many  parts  of  the  formulas 
on  the  presumption  that  they  are  too  remote  from  the  right  end  (the  end  actually 
considered)  to  be  significant.  This  latter  presunption  can  account  for  defects 
of  symmetry  about  the  origin  observable  In  some  of  the  formulas. 


* It  may  be  noted  that  limits  of  functions  fj,  ha,  kjn,  gjn,  found  will 

still  be  based  on  the  presumption  of  a stress-free  crack.  They  are  not  to  be 
confused  with  contributions  arising  simply  from  superdislocations  terminating 
on  the  crack.  Unmitigated  effects  from  such  superdls locations  could  Imply 
Infinite  stresses  on  the  crack,  but  those  will  be  missing  from  limits  of  functions 
for  which  the  surface  stresses  are  always  required  to  be  cancelled. 
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The  aim  of  the  present  work  then  Is  to  make  the  References  [A-l,A-2] 
formulas  better  founded,  more  readily  Interpretable,  and  If  possible  more 
easily  coiq>utable  by  presenting  and  examining  several  alternative  formulation 
of  them. 

Some  Auxiliary  Quantities 

It  Is  helpful  to  begin  adaptation  of  the  References  [A-l,A-2]  formulas  by 
Introducing  several  quantities  which  facilitate  their  expression  In  real  form. 

A starting  point  Is  to  observe  that  a recurring  quantity  In  their  formulas  Is 


(Zj  - a^) *,  In  which 


Zj  ■ aj  + l^e  ■ a(aj  + Aj  e ) , 


(A-4) 


where  by  definition  ^ polar  representation  of 

this  radlcAl,  one  may  seek  real  functions  f^(Xj,6)  and  f2(Xj,6)  such  that 

f j + 1 - (.j  V)  /.  . 

If  this  Is  to  hold,  then  squaring  both  sides  and  equating  their  real  parts  and 
their  Imaginary  parts  shows 

f^  - f^  “ -(1-Oj)  + ZUjXj  cos  0 + Xj  cos  20  , 

2fj^  f2  ■ ® ^ ’ 

Squaring  and  adding  these  quantities  shows 
2 2 

(fj  + f^)  - (l-Oj)  -2xj  (l-oj)(20j  cos  e + Xj  cos  20)  + j + ^ J cos  0 + xj 

which  can  be  shown  to  be  equivalent  to 

*2  ^ -2  „2 
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where 
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Dj  - {[(Xj+l-oj)  8ln0]^+  [2ajXj+(xj- l+Oj)  cos  . (A-5) 


2 2 

Rearrangement  of  the  above  * £2  also  shows 


fj-f2  - -(Xj+l-Oj)  sin  0.  8ln0+  [20jXj+(Xj  - 1 + a^)  cos  0]  cos  0, 


2 f j “-(Xj+l-Oj)  sin  0.  sin  0 + [2ajXj  + (Xj  - 1 + a^)  cos  0]  co8  0+Dj 


Let  be  the  angle  in  the  range  - ir  < 0^  £ ir  for 


which 


si.  Sj  - 


(X?  + 1 - a^)  sin  0 


2o,  X . + (X^  - 1 + a?)  cos  0 


, cos  « 


(A-6) 


Presuming  X^  ^ 0,  0 _<  Oj  _<  1 and  0 _<  0£  y so  that  sin  0^  > 0,  it  follows 


2 a,X,  + (Xj  - 1+ a^)  cos  0 


Bj  ” arccos 


(A-6’) 


It  can  be  shown  also  that  sin  (0  + 0j)  ^ 0,  so  it  Is  implied  that  0 £ 0^  £ ir  - 0. 
Now 

2 2 2 2 
2 ■ Dj  [cos  0 cos  0j  - sin  0 sin  Bj  +1]  ■ 2 Dj  cos^  — . 


2 2 2 1 
Thus  f “ + Dj  cos  , and  since  f + £2  “ Dj , it  follows  that  f2  ■ + sin  — ^ 

2 

To  have  2 f.  f.*  2a. X.  sin  0 X.  sin  20,  it  is  necessary  that  the  ambiguous  signs 
JJ  j 22^ 

for  f^  and  f^  be  alike,  and  to  have  (Zj  - a )^  in  the  right  half  plane  (as  appar- 
ently was  presumed  for  the  formulas  being  considered)  that  sign  must  be  positive. 


A-6 


Thus, 


fl  " — 


£2  - Dj  sin  -Y 


and  then 


2 2 ^ 0+6.  l^fi+R  1/7 

(zj  - a ) /a  - Dj  [cos  — + 1 sin  ~~2^  " ® 


(A-7) 


This  expression  holds  for  all  values  of  6,  A^,  and  in  the  ranges  that  were 
presumed  for  them. 

Letting  Pj^l-Uj  ■ (a-aj)a,  one  can  Interpret  by  showing 


Dj  “ » where 

2 2 ^ 

6j  “ [pj  - 2PjAj  cos  6 + Xj]  « (distance  from  Zj  to  a) /a. 


2 2 

6j  * 1(2  - Pj)  + 2X  (2-pj)  008  0+  Xj]  > distance  from  z^  to  -a) /a. 


(A-8) 


These  distances  are  Illustrated  In  Figure  A-2a . The  angle  Is  harder  to  Interpret, 

but  It  can  be  approximated  by  the  angle  between  the  lines  from  z.  to  a and  from 

Zj  to  a.  Other  quantities  of  occasional  Interest  include  and  the  polar 

coordinates  of  z.,  which  are  also  shown. 

J _7  7 


show  that 


J -2  2v** 

A second  radical  which  arises  Is  (z^  - a ) . To  reexpress  it,  one  may 


f-2  2.f  2 2.  -4  4 

(zj  - a ) (Zj  - a ) - Dj  a 


Thus, 


<ij  - -iDj  ./(.j  - -J 


(A-9) 


Since  this  rsdlcsl  too  seems  to  have  been  presumed  to  lie  In  the  right  half  plane, 
the  ssblguous  sign  here  will  be  taken  to  be  positive. 
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To  prepare  for  handling  the  quantities  (z^  - Zj)and  (z^-Zj),  let 

2 2 k 

C.  - {[(o„-o,)  + (X„- A.)  cos  e]  +[(A„-A,)  sin  6]  } , 

1 n j n j n j 

2 2 ^ 


and  let  and  ^2  angles  in  the  range  - iT<^_<irfor  which 


cos  ^ ^®n”  ^ ^^n ” 8in0  /Cj^ 

cos  “ [(a  -a,)  + (A  - A,)  cos  0]/C,,  sin  f-  “(^-  + ^4)  8ln0  /C- 
Injnj  z ^nj  z 

Then 

(z  - z )/a  - C e^*^  , (z  - i.)  - C,  e^'^2 
n j 1 n j A 

The  new  quantities  C^,  and  ^2  depend  on  both  z^  and  Zj.  They  are 

Illustrated  in  Figure  A- 2b. 


(A-10) 


(A-11) 


(A-12) 
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Functions  Related  to  Stresses 


Turn  now  to  the  functions  determined  In  References  [A-l,A-2]  to  define 
the  stress  of  five 


f ■ 4 sin  9 
n 


Re 


t 2 

(Zn-a  ) 


z - a 
n 


-fat  Re 
n 


16 

e ^ 

2 2 *S  ! 

(z  -a)(z^-a  ) 
n n 


(A-13) 


Using  results  found  above,  one  now  finds 


Re 


/ 2 2.** 

(Zn-a  ) 


z - a 
n 


- Re 


n n n n 


D -p  I 
_ °L  n 


0+B  e-B 

cos  — 5-^+ A cos — 

z n z . 


6‘ 


a 1 Re 
n 


le 


(z  -a)  (z^-a^)** 


A 

n 


,i(e-Pn)/2(-p  +A  e^S) 


n n 


r 

^nL-Pn®°®  — 


2 2 


D 6^ 
n n 


a 


Thus  the  function  f^  is  equivalent  to 


4 sin  9 

D 6^ 
n n 


[(A^  - p 
L n n n 


) cos 


e+B 


^ + A„  (D  - p ) cos  , 
2 n n n 2 


(A-14) 


This  function  represents  the  contribution  that  the  superdls locations  at  z^  and  z^ 
make  to  the  singularity  cancelling  equation.  However,  those  superdislocations 
were  presumed  to  have  senl-lnflnlte  tails  which  are  not  relevant.  Ue  can  find  an 
adjustment  to  delete  the  tall  effects  by  taking  Note  that  this  adjust- 

ment does  not  violate  boundary  conditions  on  the  crack  since  It,  like  f^.  Is 

based  on  the  presumption  that  no  stress  is  transmitted  across  the  crack.  Now  as 

2 '* 

A -K),  one  has  B ■♦p  , D_-»’l-a“  and  B •*%-%,  so 
n n n n u ii 


llm 

‘a  0 

n 


i + 0]  - 0 

p_  r l-a_  ^ 


(A-15) 


I 


A-10 


Since  this  adjustment  vanishes,  the  function  £ applies  also  when  the  superdls- 

n 

locations  are  presumed  to  lack  the  seml-lnflnlte  tall. 

The  function  f^  given  above  Is  an  abbreviated  form  based  on  the  presump- 
tion that  contributions  from  the  superdlslocatlons  at  the  left  end  of  the  crack 
are  negligible  for  the  singularity  cancelling  equation  for  the  right  end  of  the 
crack.  For  some  purposes  this  approximation  Is  imaccep table,  and  then  a fuller 
form  of  f^  based  on  an  entire  set  of  four  superdlslocatlons  Is  desirable.  To 
find  the  fuller  form  (say  f^),  one  may  use  the  form  of  the  singularity-cancelling 
equation  given  on  page  78  of  Reference  [A-1].  That  form  (for  plane  strain) 

Is 


oa 


Gb 

4v(l-v) 


4 

V 

2 2 ^ 2 2 ^ 
(z‘-a2)  (z;-a^) 

4 _ 

4 y- 

j-i  ^ 

Zj-a  ' zj-a  ^ 

j=l  ® 

, 

(A-16) 


and  appropriate  locations  z . can  be  taken  to  be  z. , z^  -z  , -z^  while  the  appropriate 

3 J-  J.*  1 J- 

i0  *“1.0  1.0  “1.0 

strengths  seem  to  be  respectively  e , -e  , e , -e  . With  these  specifications, 
one  finds 


(A-17) 


(A-18) 


As  X -»-0,  this  function  vanishes,  so  here  too  the  deletion  of  the  tall  effects 
n 

does  not  change  the  function.  It  may  be  noted  that  this  fuller  function,  though 
It  Involves  more  siq>erdlslocatlons,  has  a simpler  real  form  than  does  the  func- 
tion based  on  only  two  superdlslocatlons. 
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The  second  stress-related  function  is 


h ■ sin  0 I cos  0 + a i Re 
n I n 


„210 

e 

, 2 2.3/2 
(z  -a  ) J 


(A-19) 


f2ie  ) 

(zO)3/2  } 


^n  ^(jj3^1(3W-3fr 
n 


\ K 

i>)/2r 


so  the  real  form  of  h Is 

n 


h ■■  sin  0 I cos  0 + — cos  — — 

“ L 2 

n 


(A-20) 


This  function  represents  contributions  the  applied  load  makes  to  shear  stress 

along  the  slip  plane  at  z . Since  these  contributions  are  not  made  by  the  dls- 

n 

location,  there  are  no  tall  effects  to  be  adjusted. 

Before  considering  the  function  k , It  Is  well  to  consider  k^  for 

nn  jn 

J n.  This  Is  given  as 


V - - H 


-z^)+e“*^z  .-z^)  e^®  (z  .-z^)+e^®(z  -z^) 


• (A-21) 


Using  qtiantltles  defined  above,  one  now  finds  the  real  form 


kjn  ■ ^ cos  (0-^j^)  cos  (20-24»j^)-^  cos  (0-H'2)  (20-24>2)j 


(A-22) 


To  get  an  adjustment  for  the  undeslred  tall  effects,  note  that  as  ^^'*'0: 


~ (distance  from  a.  to  z }/a  , 

j “ 


(A-23) 


A-12 


and  t where 

1^0 


(ct  -a . )+X  cos  0 

♦«  ■ -°-'c  ° 

O 


, sin  A 


A sin  6 
n 


o C 


(A-24) 


The  quantities  and  are  Included  in  Figure  A- 2b.  Using  them,  the  adjustment 
needed  for  tall  effects  in  k.  is  seen  to  be 

Jn 


11m  . 

'\-*0  *jn 
n 


■ ■ 4 [ 


cos(e-4  )-  cos  (9+(>  ) I cos  (29-2^  ) 


4 

■ - sin  0 sin  4^  cos  (29-2^) 
o 


4A  sin  0 

-^-5 cos  (20-2(|.  ) . (A-25) 

aC^  ° 

o 


This  adjustment  for  tall  effects  in  is  not  trivial,  but  its  ultimate  effect 
is  yet  to  be  found  since  k^^  in  practice  is  always  added  to  gj^,  which  is  yet  to 
be  discussed.  We  add  the  observation  that  k.  represents  contributions  that  the 

_ J“ 

dislocations  at  and  z^  make  to  stress  on  the  slip  plane  at  z^.  if  effects  from 

dislocations  at  -Zj  and  -z^  were  included,  would  be  changed  slightly  in  its 

value  but  perhaps  more  noticeably  in  its  form. 

For  the  function  k^^,  it  has  been  shown  that 
nn 


k ■ - 7“  cos  20 
nn  t 

n 


(A-26) 


This  function  excludes  the  shear  stress  contributed  at  z^  by  the  dislocation  at 

Zn»  the  dislocation  is  not  considered  to  drive  itself  forward,  but  the 

contribution  from  the  dislocation  at  z^  is  included.  In  adjusting  for  tail 

effects,  however,  it  seems  appropriate  to  exclude  effects  from  tails  of  both 

the  dislocations  at  z^  and  z^.  We  can  find  the  adjustment  required  here  from 

that  found  for  k,  by  considering  there  a case  for  which  a,  ■ a and  I >0  but 

j n n 

letting  £.-*0.  Since  then  C -♦•A  and  4 ■♦0,  it  follow  that 
j o n o 


tall  adjustment  for  k 


_ _ 4 sin  0 


(A-27) 


Again  the  ultimate  effect  of  this  adjustment  is  tied  to  an  adjustment  to  g which 

no 

is  yet  to  be  diecueeeds 
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Before  considering  the  function  g , It  Is  well  to  consider  the 

nn 

function  gj^  for  J n.  This  is 

g - 2 sin  Re  — 5 — 5-r  — ^ 5 

f 

[,  2 2.h  ,-2  ^2.Ji-l 

v»j  V'j  J 

n 

/ 2 2.Js  L .2,  2 2.%  2 2 ^ z-z,  J 

(z^-a  ) I (z^-Zj)  (Zj”a  ) L “ J 


(A-28) 


2 2k  ,-2  2.%^  2 
-a  r (z.-a  ) + a - z,z_ 


® ^ S_J.  +2  — ^ 

T - x2,-2  2k  2 2 z-z, 

(z  -Zj)  (Zj-a  ) L z -a  n j 


Using  quantities  defined  above,  this  g^^  can  be  expressed  in  the  real  form: 


2X  sln^e  D,  40+e,-B  -4^.  D D 20-6  -B  -4(()  D 

g -—5— -^coB i.-f  cos  (26-2(1.  ) +4 cos ^ jcos  (20-2(^  ) 

*3“  ® °n  ^ cj  ^ 2 ^2  2 


a D,  20f6  -36  -2<^,  a D,  6 +36  +2(|.  X D 40+6.-36  -2(|.. 

+ -Mco. i °-  j^.f-S4coa  d > ° ^ 
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4e-e,-3  -6(f>  60-3,-3  -64., 

q)  ^~2 V ‘'°®  2 


86-3,-3  - 


n~^*l 


X,  3,-33-24,  2X  D 

— -^—5-  cos  ^ + — cos  (6-34,) 

Wn  ‘=2 


(A-29) 


[(1-0,0  -X  X 
_ j n j n 


20+3,-3  -64,  46+3,-3  -64, 

) cos— 4S_^-a^X„co. c.s-l-|-^_ 


To  find  the  adjustment  needed  to  cancel  effects  from  the  unwanted  tails  of  the 
dislocations  with  branch  points  at  and  Zj,  note  that  If  then  C,^->C2'KJq, 

Dj-^jq“(1-Oj)**  and  3j-*^-0*  Then  the  desired  adjustment  proves  to  be 


4X  sin  0 


■o  *jn 


cos (26-24^) 


(A-30) 


This  adjustment  too  Is  not  trivial,  but  It  Is  the  negative  of  the  adjxistment  for 
kj^,  so  the  adjustment  needed  for  + 1b  zero.  Since  g^^  and  are  always 
summed  In  being  used,  this  implies  that  the  unadjusted  formulas  for  them  can  be 
used  In  practice,  even  though  they  are  questionable  when  viewed  Individually. 

It  should  be  observed  also  that  gj^  as  shown  here  refers  only  to  contributions  to 
slip  plane  shear  stresses  at  z^  arising  from  forces  required  on  the  crack  to 
abolish  the  stresses  caused  by  the  dislocations  at  Zj  and  Zj.  If  effects  of  the 
dislocations  at  -Zj  and  -z^  were  Included,  there  would  be  small  additions  to  *Jn- 
In  view  of  the  long  real  expression  for  g^^,  it  is  to  be  hoped  that  such  addition 
would  provide  a 8liq>ler  overall  sum,  as  may  well  be  possible. 

Finally,  turn  to  g^^.  This  Is 


g >21  sin  6 Re 
nn  n 


r ( 2.;^. 
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« 2ie  / 2 le  ..  r 2,  . . , 2 2-Ss  ,_2  2>s^  2 _ t / 

i e ( a z e -16  fa  (z  -z  ) (z  -a  ) (z  -a  ) +a  -z  z Tj 

a / n ■ e nn.^^n  '‘n  n n 1 

-2  2.k)~  , 2 2.512  ~ 7~J~2~J^\  2 2^  z-z  ( 

(z  -a  ) I (z  -a  ) (z  -z  ) (z  -a  ) L z.~a  n n J ) 


n n n 


(A-31) 


Again  using  quantities  defined  above,  this  g^^  can  be  expressed  in  the  real  form 


2A  sin  ^0  r (2a^+l)  2a  A A^  o sin  26 

V”  26„ |-^co8  (e-26„)  --ycos  (2e-2B„)  - ^ ° 

d2  2d^  n p2  n jj2  n 2A^  sin  0 

n •-  n n n 


nn 


sin  ( 6-26  ) a A 

1 XX  XX  XX 

+ — s — : — :: — 7—  cos  36 


2 sin  0 


n 


A^  sin  6 D^  + l-o^l 

1 - -f  cos  (0-36^) 2“^  - -a ^ 

“ " 2D^sin0  2A^  J 


(A-32) 


One  cannot  manipulate  this  function  Itself  to  find  Its  unwanted  tall  effects,  but 

an  adjustment  for  them  can  be  found  by  beginning  with  g.  and  assuming  a.  > a , 

ja  J n 

but  taking  A^  > Aj.  This  assumption  alone  makes  “ 6»  aad 

thereby  simplifies  g,  . Continuing  by  letting  A, ->-0,  so  that  also  C ->C,-»-A  , 

ja  L j 1 ^ n 

D.-H),  ■ (1  - a2)^  and  6.-Hr-6,  one  finds  the  adjustment 

2 j jo  n j 


11m 
Aj-*-0  *jn  " 


4 sln^  6 


“ adjustment  for  g^^. 


(A-33) 


This  adjustment  Is  the  negative  of  that  required  for  k 


nn 


Since  g and  k are 
nn  nn 


always  added  In  use,  their  adjustments  cancel  each  other.  However,  here  too  the 

effects  of  the  dislocations  at  -z  and  -z  were  omitted,  and  Including  them  coulc 

n n 

make  minor  changes  In  value  and  perhaps  substantial  changes  In  form. 


Discussion  of  the  Functions  Belated  to  Stresses 

Up  to  this  point,  the  effort  to  make  the  stress  related  functions 

<^n'  “n*  ^n*  '^nn*  «jn*  « more  manageable  has  provided  real  expressions  for 

them  and  has  shown  how  they  would  be  changed  by  elimination  of  effects  from 
the  unwanted  tails  of  the  dislocations.  The  real  expressions  are  given  In 


A-16 


Equations  (14),  (20),  (22),  (26),  (29),  and  (32),  respectively,  using  new 

symbols  defined  In  Equations  (4)  to  (8),  (10),  (11),  (23),  and  (24) 

and  Interpreted  In  Figures  1 and  2.  Nontrivial  adjustments  for  tall  effects 


were  found  for  the  last  four  functions,  as  shown  by  Equations  (25),  (27),  (30), 
and  (33),  but  these  adjustments  cancel  each  other  when  and  gj^  (or  k^^  and 
g^^)  are  added,  as  they  always  are  In  actual  us**.  It  was  observed  that  the 
given  stress  related  functions  are  based  on  ^ .ots  from  the  pair  of  dislocations 
at  one  end  of  the  crack,  but  not  from  the  pair  at  the  other  end  since  they  were 
regarded  as  remote.  For  one  function  (^q)»  revised  forms  were  deduced  to  Include 
effects  from  all  four  dislocations  of  the  s3nBBetry  set.  Surprisingly,  the 
form  of  f^  based  on  four  dislocations  [as  In  Equations  (17)  and  (18)]  Is  simpler 
than  that  based  on  two  [as  In  Equations  (13)  and  (14)],  though  that  might  not  be 
true  for  the  other  functions. 

The  complexity  of  these  stress  related  functions  prompts  a desire  for 
sliiq>llfled  forms,  especially  for  g^^,  but  care  must  be  exercised  that  any  sliiq>llfled 
form  Is  at  least  approximately  valid  over  the  range  of  Its  parameters  for  which  It 


Is  used.  The  functions  differ  In  complexity,  so  they  will  be  considered  separately. 

The  aim,  of  course,  can  be  either  to  make  the  formulas  more  comprehensible  or 

to  make  them  more  rapidly  computable. 

The  form  of  f based  on  four  dislocations  [as  In  Equation  (18)]  Is  the 
n 

sliiq>ler  one  for  f^,  so  consider  it  first.  It  has  three  terms,  of  which  the  first 

two  have  essentially  the  same  magnitude  while  the  third  Is  smaller,  since  and 

X /D^  have  the  order  of  unity  while  X is  smaller.  To  combine  the  terms  of  highest 
n n n 

order,  we  may  use  the  Identity 


A-17 


1-a  1-0 

Here  the  first  term  should  be  dominant,  since  X and  — ) are  expected 

XQ  A 

n n n 

to  be  small,  while  — r has  the  order  of  unity.  Indeed  If  (l-o  )/X  Is  small 
2D^  n n 

(that  Is  (a-a^)/!^  Is^small),  then  from  Equation  (5) 


D»/2X  , 

n n 


and  from  equation  (6) 


X 1-0 

®n  ^~2  ^ ~x~°^  ® ’ 

n 


which  Is  small.  It  can  be  seen  that  the  first  order  approximation  Is 


fT  , fl  . o 0 

f w 6 /t—  sin  0 cos  — , 

n A 2 

n 


(A-36) 


In  agreement  with  the  approximation  given  In  References  [A-l,A-2].  If,  however, 

one  wishes  to  account  somewhat  for  the  distance  l-o  (=p  ),  some  of  Its  effect 

n n 

can  be  seen  In  each  of  the  terms  In  the  brackets  of  Equation  (35)  as  well  as  In 

D . Letting 
n 

e “ P„/^_  “ * 

n n n n n 

and  keeping  terms  of  order  X^  and  e^,  one  finds  from  Equations  (5)  and  (6) : 


X +2e 

B « - ° sin  0 . 

n 2 
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Using  these  approximations.  Equation  (35)  yields  the  approximation 


sin  6 cos 


A 5e 

if  - (1-2  cos  0) 


(A-37) 


T 

This  approximation  to  provides  some  insight  concerning  the  relative  effect  of 

the  distance  p under  the  presumption  that  A and  p /A  are  small,  though  its 
u n n n 

computational  advantage  over  the  exact  form  in  Equation  (35)  may  be  small. 

It  is  Important  now  to  note  that  if  the  function  f as  given  by  Equation 

T ^ 

(14)  is  approximated  to  the  same  degree  as  is  in  Equation  (37),  one  finds 


fjj  = e/^  sin  e cos  I [l+^  - (1-2  cos  0)]  . (A-38) 

T 

This  approximation  differs  from  that  for  f by  a term  of  order  A In  the  brackets* 

n n 

The  difference  between  these  approximations  must  be  due  to  the  difference  between 

T 

the  original  complex  representations  for  f and  f as  given  by  Equations  (13) 

n n 

and  (17).  Since  Equation  (17)  has  better  jtistlficatlon  (from  using  all  four 
dislocations  of  the  S3nmaetry  set  Instead  of  only  two) , the  approximation  (37) 
must  be  regarded  as  better  than  (37).  Thus,  :*n  this  sense.  Equation  (13)  does 
not  support  the  degree  of  approximation  shown  in  Equation  (38).  Since  the  com- 
plex representations  of  g^^,  and  g^  as  given  originally  by  equations 

(21),  (26),  (28),  and  (31)  are  based  on  the  same  neglect  of  the  two  more  remote 
dislocations,  it  must  be  considered  doubtful  that  they  would  support  a higher 
order  approximation  as  typified  by  Equation  (37),  though  they  should  support 
first  order  approximations  as  typified  by  Equation  (36).  Since  the  first  Influence 
of  is  the  same  in  both  Equations  (37)  and  (38),  it  might  be  hoped  that  the 

first  influence  of  e could  be  correct  for  the  other  stress  related  functions 

n 

even  thou^  the  formulas  are  based  on  two  dislocations,  but  use  of  four-dis- 
location forms  would  make  calculations  more  secure. 

To  get  an  approximate  form  of  h^  from  Equation  (20)  to  accuracy  comparable 

to  that  in  Eqtiatlon  (37),  one  may  reuse  the  same  approximations  for  D and  B 

n n 

and  thus  show 
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sin  0 cos  0 + cos  | [l  + (1-2  cos  0)  + . (A- 39) 

The  first  two  terms  of  this  approximation  agree  with  those  given  in  References 
[A-l,A-2].  Since  this  function  represents  effects  from  the  applied  load  which 
are  not  dependent  on  how  many  dislocations  were  employed,  the  higher  order 
terms  here  should  also  be  valid. 

Approximation  of  the  function  given  by  Equation  (22)  involves  con- 
sideration of  more  parameters  since  there  are  two  sets  of  dislocations.  In  order 
to  express  the  results  more  compactly,  let 


A,  - a - o , L . ■ X - 1..  L “X  + X, 
dnj’dnj*snj 


l/slng  these  quantities,  one  may  express  the  component  parts  of  Equation  (22) 
exactly  and  thus  rewrite  that  equation  in  the  form 


2 2 
sin  e 


— ^ [A^  cos  8 + L^cos^  ® ” ^8  ® ^ ^d  ® 


e]^-i| 


where 


2 2 2 2 2 -2 
C,  - A^  -f  2 A ,L  . cos  0 + Lj  cos  0 + L sin 
2 d d d d s 
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/ 

I 


To  simplify  this  result.  It  will  be  assumed  that  A.  Is  of  smaller  order  than 

d 

or  Lj.  Ej^andlng  the  various  parts  of  the  expression  for  as  series  In 
(Aj/Lg)  or  truncated  after  terms  of  order  1,  one  finds  at  last  that 


2 2(L^  cos^6  - Lg  sin  ^6)  [2(L^  cos^0  + sln^  ^s  ^ 

“’^Jn  ■ IT  .,2  2 . . ,2  . 


(Lj  cos^  0 + sin  ^0) 
d 8 


- 2A , cos  0 
a 


- 2a,co8^0+L  8ln^0)(L.  cos^O-  3L  sln^^+dj+L  :?8ln^0-L5 

1 d 8 ''d  s ds'  d 


Ur 


/t2  2a  ^ ,2  . 2„v2 

(L . cos  0 + L sin  0) 
d 8 


] 


fA-AO'i 


It  can  be  seen  that  If  X . 0 so  that  L . ■ L “X  , this  equation  reduces  to 


L » n J X 

n 

a form  In  which  the  term  -cos  20/X^  Is  like  one  given  In  References  [A— l,A-2], 

while  the  term  1/X  reminds  one  of  the  deleted  self-driving  effects  In  k . 

n . nn 

There  was  no  term  corresponding  to  that  Involving  A,.  Thus,  the  k^  of  Equation 

^ jn 

(40)  Is  fairly  plausible.  The  accuracy  of  approximations  using  the  second 

order  term  may  be  doubtful,  however,  since  the  two  more  remote  dislocations  were 

neglected  In  the  original  formula  (21),  though  the  contribution  of  the  term 

Involving  A^  may  be  nearly  right. 

The  value  of  k given  by  Equation  (26)  needs  no  approximation, 
nn 

The  function  of  g^^  as  expressed  by  Equation  (29)  contains  all  the 

quantities  appearing  In  the  preceding  functions.  Including  two  distances  D and 

angles  0 corresponding  to  the  two  points  and  z^.  Some  formal  sliiq>llflcatlon8 

of  the  function  can  be  had  by  Introducing  also  (R^,  ^j)  and  (R^,  ^q)>  the  polar 

coordinates  of  z.  and  z . Thus,  g.  can  be  written  as 
3 n ]n 


2X  sln^  0 

- . _J! 

*jn  aD 


cos 

"l 


40+8.-6-40, 

i-j-S 008(20-24^)  + 

^1 


+ -A  cos 


20,8^-8^-4»2 


“2  cos  (26  - 24«)  + 
I 


(A-41) 
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D-R  20+e  -3e„-2^. +2* 
in  .1  n X n 

-L-J  COS 2 ^ 

C,D‘‘ 

1 n 


Va  ' 


A 26-0  -36  -2(>  2A  D 

^ cos  J-5-S i + — La  cos  (3e-3«|.,) 

C,DjDJ  <=1 


2v  r 

+ COS 

^pj 


4e-s-6„-6*i 


- R R.  cos 
n j 


4e-6£^-6(|»i  + 2t|>  +2i|) 


A,  6,-36  -2*-  2A  D 

— -J— r cos  ^ \ + ~ 4 ° cos  (e-S^,) 


r 

+ cos 

L 


2e+6.-6jj"^®2 


- R R,  cos 

n J 


26+6,-6  + 2t|>  - 2* 

1 n X n 1 


It  Is  helpful  to  note  that  all  the  angles  6j,  6^.  ♦j,  and  are  small  to  the 

order  of  the  relevant  A and  e(>  p/A).  The  tvo  bracketed  quantities  In  Equation 

(41)  are  also  small  to  this  order,  so  that  all  terms  In  the  braces  have  order 
-3/2 

approximately  like  A 

Approximation  of  the  function  for  small  values  of  selected  para- 
meters Is  much  more  tedious  than  the  process  for  the  earlier  stress-related 
functions,  since  g^^  Involves  so  many  parameters,  and  because  It  has  so  many 
terms.  Including  a great  variety  of  trlgommetrlc  functions.  In  order  to  make 
the  approximation  possible,  all  these  parameters  were  treated  as  being  small 
to  the  same  order:  Aj,  A^,  L^,  L^,  6jt  3^^*  ♦j*  ♦jj*  ^d^^d 

treated  as  small  to  higher  order).  In  addition,  a new  quantity  was  defined  to  be 


- (lJ  cos^  e + sln^  e)** 

c d 9 


so  that  It  also  Is  small  to  the  same  order.  With  these  assumptions  and  much 
algebra,  it  was  found  that  the  first  approximation  to  gj^  is 
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that  X.  » X or  that  X » X..  Here,  by  taking  a.  ■ a (so  that  A.  ■ 0)  and 
j n ^ j j ^ ® 

Xj  > X^,  one  finds  from  Equations  (40)  and  (42)  that 


2(1-  cos  46) 


6 

2 


(A-45) 


Adding  these  results  provides  the  same  approximation  as  Is  shown  for  ^X2  °° 

Page  37  of  Reference  [A-1].  The  case  with  a.  ~ a and  X > X Is  somewhat  more 

j n n j 

delicate,  since  the  leading  terms  of  the  respective  approximations  cancel  each 

other.  Therefore,  It  Is  necessary  to  keep  terms  of  the  first  two  orders  In  X^. 

[This  Is  not  to  be  confused  with  the  first  two  orders  over  all  the  parameters, 

as  governed  the  approximations  leading  to  Equations  (4)  and  (42).  This  Is  merely 

a way  of  selecting  the  remaining  major  terms  In  those  two  eqtiatlons.]  With  this  In 

view.  Equations  (40)  and  (42)  yield  for  X » X. . 

n 1 


*111 


n 


*Jn 


2 

• - ^ ^ ^ ^ (l  + cose)sln^  0 

n 21 


(A-46) 


In  the  sum  k.  g.  , the  leading  terms  here  cancel  each  other,  and  the  resulting 
Jn  jn 

sum  Is  the  same  approximation  as  Is  shown  for  k2j^  + g2]^  on  Page  37  of  Reference 
[A-1].  The  leading  terms  which  here  cancel  each  other  are,  of  course,  the 
quantities  which  would  have  been  removed  alternatively  by  the  coiq>ensatlon  for 
tall  effects  In  the  fmctlons  k^^  and  gj^. 

It  may  be  observed  also  that  the  distinction  between  the  approximations 
offered  by  Equations  (40)  and  (42)  and  those  cited  from  Reference  [A-1]  consists 
mainly  In  the  more  general  lengths  1^  and  1^  allowed  by  Equations  (40)  and  (42). 


It  is  Inpliclt  here  that  accounting  for  finite  values  of  both  these  lengths  is 

more  Important  than  is  the  separation  between  the  slip  planes,  which  determines 

Aj”  (a  -a,) /a.  This  order  of  Importance  is  quite  plausible  except  if  i.  and  I 
3 i Q 

are  nearly  equal,  but  such  cases  were  ruled  out  for  Equations  (AO)  and  (42)  when 

it  was  assumed  that  A,/L , is  also  small,  that  is  that  (a  is  small. 

o G n j n j 

This  assumption  should  be  remembered  in  applying  these  approximations. 

The  great  abundance  of  parameters  entering  these  stress-related  fimctlons 
places  a great  burden  on  the  organization  of  any  approximation  process.  This  is 
well  Illustrated  by  the  complexity  of  the  second  order  terms  of  g.  given  later. 
Some  simplification  might  yet  be  found  for  the  sum  of  those  terms,  but  it  probably 

would  be  minor  unless  some  other  drastic  assumption  is  made,  such  as 
taking  Xj  » X^  or  X^  » X^ . Assump  tions  such  as  those  could  be  imposed  readily 
on  the  second  order  terms  for  g^^,  but  they  probably  would  still  yield  terms  com- 
peting in  magnitude  with  the  discrepancy  allowed  when  the  two  remote  dislocations 
were  discounted.  This  Illustrates  the  caution  needed  in  approximation  of  these 
functions.  It  is  good  that  exact,  real  expressions  exist  for  all  the  stress-related 
functions,  and  but  for  gj^  they  are  siiiq>le. 

Formulas  for  Crack-Opening  Displacements 

The  Equation  (3)  for  the  crack-opening  displacement  employs  a function 
dj  for  which  Reference  [A-1]  gives  the  complex  form 


J ( ra^-xz  +I(z5-a^)**(a^-x^)‘*-j  ra^+x2 

* » ^ r°‘  L.^-«  ^ 


r.2.  2 2.Js-  2 2.% 

^ +x^  + i(z~-a  ) (a  -X  ) 


(A-47) 


This  d.  multiplied  by  b.  sin  6 represents  the  contribution  to  the  crack-opening 
J ^22 

displacement  at  position  x (for  x ^ a ) made  by  the  set  of  four  dislocations 
symmetric  to  the  one  with  branch  point  at  z^,  each  extending  semi-inf Initely. 
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To  begin  to  find  In  real  form,  one  may  show  first  that  In  terms  of 
the  qiiantitles  Introduced  above  and  the  position  variable 

one  finds 


I lj(a^-x^)*Re 


¥ 


le 


>2  2 55  , 2 2. 

(Zj-a  ) (Zj-x  ) 


(A-49) 


e+3B.  e-3e 

cos  — *** 


] 


0-e.  30-g^ 


cos- 


vd: 


[1  + 2 vj  cos  (0  + 8j)  + Vj] 


+ cos 


To  prepare  for  treating  the  logarithmic  terms  In  dj,  consider  any  real  numbers 
a and  b,  and  define  ^ (where  - ir  <♦£  ir)  so  that 


cos  ^ ■ 


/Xl2 

/ a +b 


sin  ^ 


/J^ 


that  Is  so  that 


Then 


^ • sgn  (b)  arccos 


/ 2^2 
/ a +b 


and 


log  (a-Kb)  ■ log  (cos  t + 1 sin  + log 

“ log  e^®  + log  1^^+b^  ■ 1 ^ + log  *4^4b^ 


Im  [log(a*Mb)]  ■ ♦ ■ sgn  (b)  arccos  — * 

4v 


(A-50) 
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Applying  this  to  terms  to  appear  In  d^,  and  letting  s " x/a,  one  finds 


I ra^-x2,+l(z^-a^)**  (a^-x^)**1  ra^+xz.+lCz^-a^)**  (a^-x^)**  I | 

^ ^ \ / 2 2 % , 2 2.^8  l''  2 , 2 2A/  2 2Ji  I 

( La  -xZj-l(Zj-a  ) (a  -X  ) J \_a  +xZj-l(Zj-a  p (a  -x  r J ) 


a^+l(zj-a^)?*(a^-x^)fs 

‘”1}  1 

■ 2_  - 2 2.J5  - 2 2Js’ 

a +i(z.-a  (a  -x  Y 

+ xzj 

V-l(z2-a2)’*(a2.x2)’s- 

•”1}  i 

■ 2 . , 2 . 2 J5  , 2 2^ 

a -i(Zj+a  >* (a  -x  Y 

i^-H(z^-a^)^  (a^-x^)^]  ^-x^zj  | ^ I (a^-x^)+21  (a^-x^)^(z^-a^)^+l^  (z^-a^)  | 


ir  I f 2 . . 2 2.*S  . 2 2v%  12  2 2 1 ir  o 2.  . 2 2 58.  2 2A..2,  2 2J 

( j^a  -l(Zj-a  ) (a  -X  ) J -X  z^  ) j (a  -x  )-21(a  -x  ) (Zj-a  ) +1  (a^-a  )) 


~ Im  J log 


— Im  j log 


r (a  -X  ) +l(z.-a^rl  i ^ I r (a  -X  ) + 1 a(f  j^+1  £2)  "I 

\ , 2 ZM  ,,  2 2 ’sj  I*  I®  j log  2 2 >S  . f X 

L(a  -X  ) -l(Zj~a  ),-J  ) j L (a  -x  ) -la(fj^+lf2)  J 

(*C?-f2)+lf^  I 2 r /LV-f2 

— , y I " "ir  arccos  + arccos  - 

(/i-s^'+f  )-if  r~ri.  2 2 

^ I L J /yi.-  _ f \ J.  <:*  J 


/ 1-s^-f 


! T ^2  / / 2 2 2 

(/l-s  -£2)  +£^  / (/1-s^+£2)  +£‘ 


— arccos 


^(l-s^-£j-£^)^+4£j(l-s^) 


— arccos 

Tt 


, 2 ,.2.,  2 0+e 

'(v  -1)  +4v  cos  — ^ 


— arccos 

7l 


f4  2 

V V2v^  cos  (0+0  )+l 


, where  v ■ / l-s^ ^ 


(A-51) 


2 2 

For  all  relevant  a,  6,  and  x (that  Is  x £ a ),  this  last  expression  Is  a 
continuous  function  which  should  be  valid  and  which  avoids  the  many-valued 
nature  of  the  original  logarltham.  This  Is  the  function  Intended  In  the  Reference 
[A-1]  analysis.  The  variable  v appears  to  be  a natural  variable  to  replace  x, 
and  it  can  be  seen  that  it  varies  nearly  in  proportion  to  (a-x)^^2  ^ 
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Combining  the  above  results,  one  finds  that  Equation  (47)  is 
equivalent  to 


d - 26(a  -jx])  - 
■’  irD. 


41  ^Vj  Oj  cos 


0+36  0-36, 

— -j-^+^cos — 2" 


11  2r 

+ v^  g^coe — 2"*'  — 2' 


[l+2VjCOS  (0+6j)+Vj] 


1-v, 


(A-52) 


— arccos 


/1+2  Vj  cos  (e+Pj)  +Vj 


To  find  adjiistment  to  compensate  for  the  unwanted  tails  of  the  dislocations,  we 

2 2 

note  that  as  6j-*-v-0,  and 


here  H(aj  - |x|)  is  the  Heaviside  step  fimctlon.  Adding  this  adjustment  to  dj 
gives  the  conpensated  form 


4X  V g.  cos  2 . 

**10  “ 2H(a  -|x|) ^ 2 4 

J xDj  [l  + 2Vj  cos  (0+  6j)  + vp 


0+36,  0-  36. -1  2 r 

-^+A,  cos — Y’  1° 


0-6, 


30-6, 


g^  cos  2 CO®  — 2' 


1- V 


— arccos 

IT 


A. 


/ l + 2VjCOs  (0  + 6j)+Vj 


(A-53) 


To  this,  we  may  add  the  observation  that  if  v^  < 1,  then 


1-v, 


- — arccos 

ll 


/l  + 2vjco8  (0+6j)+Vj 


arctan 


0+6, 


2 v^  cos  — y' 
1-V? 


(A-54) 


A-28 


This  latter  form  may  have  coii9>utatlonal  advantages,  but  It  has  an  unwanted  step 


If  Vj  Increases  through  1. 

The  desired  ftmctlon  d^  Is  represented  exactly  by  Equation  (53),  but 

jc 


if  one  wishes  to  examine  the  behavior  of  d^  for  small  values  of  v^,  say 


|vj|<|  (that  Is  |x|>aj),  he  may  alter  dj  by  using  the  expansions 


cos  a*+v  cos  g* 


1 + 2v^co8  (a'  + 6’)+v^  n»o 


- I (-v^)  cos  [(n  + 1)  o'  + B'l  • and 


2 V cos 
- 1-v^  - 


arctan 


(A-55) 


2 V Z 


n*o 


(-V  ) -08 

2n  + l 2 


From  these.  It  follows 


« 2 n r 

J TTD  n»o  L J 


(2iH-l)0+(2n+3)3 


cos 


iTDj  n“o 


i + Xj  cos 


(2n-l)e+(2n+3)B 


] 


Av,  » (-vj)"  (2n+l)(0+B  ) 

, ' -dr  “• — 2 — ^ 

n»o 


, for  |vjl<l 


(A-56) 


To  get  simple  approximation  for  d^,  take  the  first  two  terms  of  each  of  the  series. 


This  shows 


4v 


djc-2H(aj-|x|)-  — 


|[ 


0+B.  a.\.  0+3B.  Aj 


e-3B, 


cos  — S-*-  + —S**-  cos  — =— *■  + -d  cos 


(A-57) 


cos 


36+33,  o.X 
— cos 


36+53,  xf  6+3b, 

I I.  J ^ ^ COfi  ” M 

2 2 2 


J 


j..,. 
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Inspection  shows  that  the  first  bracketed  term  here  Is  equal  to  f (8  sin  9) . 

This  Is  precisely  the  relationship  needed  to  make  the  terms  of  order  cancel 

In  the  equation  for  &u.  provided  a and  b are  related  by  the  singularity  cancelling 

^ T 

equation.  Satisfaction  of  this  relation  shows  again  that  f^  Is  a more  satisfactory 
function  than  f^.  To  arrange  the  bracketed  quantities  In  Equation  (57)  Into  terms 
of  descending  magnitude,  use  Eqxiatlon  (34)  and  the  relation 


a.  A.  3&4-58 

cos 


^ . 3W-3B.  aJ  9+58.  1-a? 

2 ^ ^ cos  ^ cos  — ^ cos 


59+58 


i 


(A-58) 


Then  equation  (57)  becomes 


4v 


dje  - 2H(aj-lx|)-  ^ 


3 9+8.  A^  9-38.  1-a^  39+36. 

^ cos  — ^ cos  — r— ^ ^ cos  

^ 2Dj  ^ 2Dj 


- V, 


6 

39+38^ 

A^ 

+ ^ 

39+58 

l-g^ 

—4  cos 

59+58^ 

2 

+ ^ cos 

2D. 

2 + 

2D. 

2 

^ ♦ 

j 

j 

Then,  assuming  that  Aj  and  Cj  (Hp^/A^)  are  small,  an  approximation  Including 
effects  up  to  levels  Included  In  Equation  (37)  Is 

djc  * 2H(aj-  |x|) ^ cos  j (1-3  cos  9)  - (5  - 7 cos  9)j  (A-59) 


- 4 (2  cos  9-1)  +-^  (2  cos  9 

Jo  « 


2 cos  78-2  cos  6-Hl 


prwwi4mt  cbeC  also  v^  « 1. 


Ki)  r 2-5  ^ °i 

L 3- A sin  9 J 


5 1-  cos  3 9*1/  1 , 

-:r  (2  cos  9-1) 5-  > ... 

^ 3-4sln^9j)  J 
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A particular  value  of  d.  which  may  be  of  Interest  Is  that  for  which 
+ ^ 

X a^,  this  being  the  position  for  which  dj  should  be  smallest.  Assuming  that 
and  Cj  are  small,  we  see  that 


Dj  = 2 Xj  and  l-a^  = 2 (1  - ) - 2 X^  * , so  that 


& 


+ 2 

Letting  X -►  a makes  v.  Using  this  value  in  Equation  (55)  shows 

j j j 


llm 
+ d. 

x-^aj  jc 


^ 0 TS  ^1 

1 n - •»  rn«  G’V J- 

2 L2 


cos 


2 (1-3  cos  0)-  ^2  ” cos0)j  = 


6 e; 


— ^ cos  -r-  . (A-60) 
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This  limit,  multiplied  by  b^  sin  0,  represents  the  Interference  between  the 
surfaces  which  this  dislocation  would  cause  In  the  absence  of  an  applied  load 
and  with  no  other  dislocations. 

The  function  d.  Is  the  one  that  Is  needed  to  define  Equation  (3) 

J ^ 

which  Is  used  to  represent  crack-opening  displacement.  We  have  here  provided 
an  exact,  real  representation  of  It  In  Equation  (53),  and  this  representation 
is  one  which  should  be  readily  coiq>utable  as  well  as  moderately  well  coiiq>re- 
henslble.  An  alternative  series  representation  of  the  exact  form  Is  given  by 
Equation  (56),  and  It  leads  to  approximations  such  as  that  offered  by  Equation 
(59).  That  approximation  must  be  used  only  In  view  of  the  assu]iq>tlons  that 
Vj  and  [=(a-aj)/^j]  are  small,  so  Its  use  In  computing  must  be  restricted 
accordingly,  but  It  also  offers  opportunities  for  Interpretation  such  as  that 
In  Equation  (60). 

It  may  be  remarked  that  computations  of  d^  based  on  the  complex  Equation 
(47)  are  apt  to  produce  misplaced  steps  In  the  results,  since  the  Individual 
logarithmic  terms  can  have  steps  for  values  of  x quite  different  from  where  the 
entire  function  should  have  Its  step.  This  difficulty  has  been  overcome  In  the 
real  formulations  shown  here. 


Towards  Some  Simple  Theorems  on  the  Displacements 
Produced  on  the  Crack  Due  to  Dislocation  Crack  Interactions 


This  report  emphasizes  how  a simple  model  of  plastic  relaxation  at  a 
crack  tip  may  be  useful  In  trying  to  give  a quantitative  description  of  fatigue 
crack  growth.  It  Is  envisioned  that,  as  the  crack  progresses,  dislocations  are 
emitted  from  the  crack  tip  leaving  In  their  wake  a displacement  discontinuity 
equal  to  the  strength  of  the  super  dislocation.  That  Is,  when  dislocations 
leave  the  crack  tip,  they  leave  behind  a stress  on  the  crack.  Using  superposi- 
tion, the  resulting  displacement  can  be  added  to  the  displacement  field  produced 
by  the  opening  of  the  crack  due  to  the  applied  stress  and  residual  stress  of 
the  dislocation. 

A question  of  some  Importance  with  regard  to  the  prediction  of  a 
crack-closure  effect  In  fatigue  Is  the  effect  on  the  crack  face  displacements 
of  the  coiiq>resslve  stress  due  to  the  dislocations.  To  bound  this,  consider 
the  following  agrument.  The  compressive  stress  due  to  an  arrangement  of  dis- 
locations cannot  be  greater  than  when  the  superdis location  pairs  are  right 
at  the  crack  tip.  This  arrangement  gives 
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Solving  this  limit  problem  using  a continuous  distribution  of  dislocations,  or 
via  a direct  simple  argument,  shows  that  It  gives  a displacement  field  of 

V ■ - 2b  H [X^-x^  ] , 

so  If  |x|  < lx|,  then  v - - 2b.  Hence,  the  net  effect  of  dislocations  emanating 
from  the  crack  tip  together  with  the  applied  stress  cannot  ever  be  negative. 
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Addendum 


Sum  of  second  order  terms  for  approximation  of  [cf.  Equation  (42)]: 
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APPENDIX  B 

CCmPUTATION  OF  CRACK  OPENING  DISPLACEMENTS  WITH  AN 
INCREMENTAL  ELASTIC-PLASTIC  FINITE  ELEMENT  METHOD 

Introduction 

Crack-opening  displacements  for  a 12-lnch  wide,  0.25-lnch 
thick  panel  with  a 1.04- inch  center  crack  have  been  computed  using  an 
elastic-plastic  finite- element  computer  program.  Three  different  panel 
materials  - aluminum,  titanium  and  carbon  steel  - were  analyzed.  The 
properties  for  each  of  these  materials  are  sunmarlzed  in  Table  B-1. 

Two  types  of  stress  loads  were  applied  to  the  panels.  One 
series  of  analyses  were  made  loading  each  of  the  three  material  panels 
to  half  their  yield  strength  and  removing  the  load.  A second  series  of 
six  analyses  were  conducted  only  for  the  aluminum  material.  The  aluminum 
panels  were  cyclically  loaded  to  stress  states  of  varying  amplitude  and 
mean  stress.  The  center  crack  was  analytically  propagated  in  one  analysis 
and  held  constant  in  the  other  five  analyses.  Table  B-2  summarizes  the 
loading  conditions  for  each  of  the  six  cyclic  analyses. 

The  following  presents  a description  of  the  finite-element 
program  and  tables  of  crack-opening  displacements  for  each  analysis. 

Description  of  the  Finite-Element 
Computer  Program 

A Battelle  finite-element  computer  program  has  been  developed 
for  the  specific  purpose  of  solving  elastic-plastic  fracture  problems. 

The  code  is  an  extensively  modified  version  of  the  AXISOL  code  developed 
by  Wilson  [B^l].  The  program  uses  triangular  or  quadralateral  constant- 
strain  elements  formulated  by  the  conventional  stiffness  approach.  Each  'i 
triangular  element  is  allotted  six  degrees  of  freedom  by  the  program. 

The  quadralateral  elements  are  comprised  of  four  triangular  elements 
defined  by  the  four  comer  nodes  of  the  quadralateral  element  and  a cen- 
troidal  node.  A plane  stress  or  plana  strain  state  can  be  modeled  with 
the  program. 
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The  final  solution  for  a problem  is  obtained  by  adding  solutions 
for  discrete  load  increments.  Element  stiffnesses  are  recalculated  for 
each  load  increment.  Elastic  loading  from  an  elastic-plastic  stress  state 
can  be  modeled  with  the  program.  Linear  or  piecewise- linear  hardening  is 
currently  employed  by  the  program.  No  Bauschinger  effect  is  considered. 

The  program  can  simulate  crack  growth  by  relaxing  node  forces 
at  the  node  representing  the  crack  tip.  Five  relaxation  steps  are  em- 
ployed to  release  the  crack-tip  force.  Sie  relaxation  procedure  is  concep- 
tually that  used  by  HcClintock  [B-2]  and  by  Andersson  [B-3]. 

Computed  Crack-Ooening  Displacements 


The  following  tables  present  crack-opening  displacements  for 
the  aluminum,  titanium,  and  carbon  steel  panels  subjected  to  one  load- 
unload  cycle  and  for  aluminum  panels  subjected  to  varied  cyclic  stress 
loads.  Tables  B-3  through  B-5  list  the  displacements  for  the  aluminum, 
titanium,  and  carbon  steel  panels,  respectively.  Tables  B-6  through  B-11 
list  the  displaceiaents  for  the  cyclically  loaded  aluminum  panels. 
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TABLE  B-1. 


MATERIAL  PROPERTIES  FOR  CENTER- CRACKED  PANELS 


Material 

Yield 

Stress,  ksl 

Elastic 
Modulus,  ksl 

Plastic 
Ftodulus,  ksl 

Poisson's 

Ratio 

Aluminum 

52.5 

10,300 

309 

0.34 

Titanium 

120 

16,000 

480 

0.34 

Carbon 

Steel 

52.5 

29,500 

885 

0.30 
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TABLE  B-2.  CYCLIC  APPLIED  STRESS  LOAD  LEVELS 
FOR  ALUMINUM  CENTER- CRACKED  PANELS 


Analysis 

Number 

Mean 

Stress,  psl 

Stress 

Aa^litude,  psl 

Table  of 
Displacements 

1 

13,125 

13,125 

B-1,  B-6 

2 

19,688 

6,562 

B-7 

3 

26,250 

13,125 

B-8 

4 

19,688 

6,562 

B-9 

5 

19,688 

19,688 

B-10 

6* 

13,125 

13,125 

B-11 

* Center  crack  extended  In  this  analysis. 
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TABLE  3-3.  CRACK  OPENING  DISPLACEMENTS  FOR  A1  2024-T3  CENTER-CRACKED 

PANELS  FOR  A LOAD-UNLOAD  SEQUENCE  OBTAINED  FROM  AN  ELASTIC- 
PLASTIC  FINITE  ELEMENT  ANALYSIS  (Plane  Stress) 


a = 0.52  Inches  w = 12  Inches 


Applied 

Stress> 

Crack  Opening  Displacement  (Inches)  At  Various 
Positions  On  Crack  Face  (Inches  From  Center) 

psl 

x*0 

x-.l 

x=.2 

x=.3 

x=.4 

x=.461 

Maximum  Applied  Stress  = 

26,250  psl 

Loading  Phase 

22,875 

.00220 

.00215 

.00201 

.00176 

.00133 

.00090 

24,000 

.00232 

.00228 

.00213 

.00187 

.00142 

.00098 

25,125 

.00244 

.00240 

.00225 

.00198 

.00152 

.00107 

26,250 

.00257 

.00252 

.00237 

.00209 

.00161 

.00115 

Unloading  Phase 

22,495 

.00221 

.00217 

.00204 

.00180 

.00140 

.00101 

18,746 

.00186 

.00182 

.00171 

.00151 

.00118 

.00086 

14,997 

.00150 

.00147 

.00138 

.00123 

.00097 

.00072 

11,248 

.00114 

.00112 

.00106 

.00094 

.00075 

.00058 

7,499 

.00078 

.00077 

.00073 

.00065 

.00054 

.00043 

3,749 

.00042 

.00041 

.00040 

.00037 

.00032 

.00029 

0 

.00006 

.00007 

.00007 

.00008 

.00011 

.00015 
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TABLE  B-4.  CRACK  OPENING  DISPLACEMENTS  FOR  T1-6AJI-4V  CENTER- CRACKED 

PANELS  FOR  A LQAD-UNLOAD  SEQUENCE  OBTAINED  FROM  AN  ELASTIC- 
PUSTIC  FINITE  ELEMENT  ANALYSIS  (Plane  Stress) 


a - 0.52  Inches  w = 12  Inches 

Applied  Crack  Opening  Displacement  (Inches)  At  Various 

Stress,  Positions  On  Crack  Face  (Inches  from  Center) 

psl 

x=,2  x*.3  x®.4  x“.461 


Maximum  Applied  Stress  = 60,000  psl 


Loading  Phase 


52,417 

.00325 

.00318 

.00298 

.00261 

.00197 

.00134 

54,945 

.00343 

.00336 

.00315 

.00276 

.00211 

.00146 

57,472 

.00362 

.00354 

.00332 

.00292 

.00224 

.00158 

60,000 

.00380 

.00372 

.00350 

.00308 

.00238 

.00170 

Unloading  Phase 


51,424 

.00327 

.00321 

.00301 

.00266 

.00206 

.00149 

42,854 

.00274 

.00269 

.00253 

.00223 

.00175 

.00128 

34,283 

.00221 

.00217 

.00204 

.00181 

.00143 

.00106 

25,712 

.00168 

.00165 

.00156 

.00139 

.00111 

.00085 

17,141 

.00115 

.00113 

.00107 

.00096 

.00079 

.00064 

8,571 

.00062 

.00061 

.00059 

.00054 

.00048 

.00043 

0 

.00009 

.00010 

.00010 

.00012 

.00016 

.00021 
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TABLE  B-5.  CRACK  OPENING  DISPLACEMENTS  FOR  CARBON  STEEL  CENTER- CRACKED 
PANELS  FOR  A LGAD-UNLOAD  SEQUENCE  OBTAINED  FROM  AM  ELASTIC- 
PIASTIC  FINITE  ELEMENT  AN/^YSIS  (Plane  Stress) 


Applied 

Stress, 

psi 

a ■ 

0.52  Inches  w ■ 12  Inches 

Crack  Opening  Displacement  (Inches)  At  Various 
Positions  On  Crack  Face  (Inthes  from  Center) 

x-0 

3^.1 

X-.2 

X-.3 

3^.4 

3C“.461 

Maximum  Applied  Stress  * 

26,250  psl 

Loadlna  Phase 

18,712 

.00045 

.00044 

.00041 

.00036 

.00026 

.00018 

21,225 

.00050 

.00049 

.00045 

.00039 

.00029 

.00019 

23,737 

.00055 

.00054 

.00051 

.00044 

.00032 

.00020 

26,250 

.00061 

.00059 

.00055 

.00048 

.00035 

.00022 

Unloadine  Phase 

22,496 

.00051 

.00050 

.00047 

.00040 

.00029 

.00019 

18,746 

.00042 

.00041 

.00039 

.00033 

.00023 

.00015 

14,997 

.00033 

.00032 

.00030 

.00025 

.00018 

.00011 

11,247 

.00024 

.00023 

.00022 

.00018 

.00013 

.00007 

7,499 

.00014 

.00014 

.00013 

.00011 

.00007 

.00003 

3,749 

.00052 

.00052 

.00005 

.00003 

.00002 

.00001 

0 

.00004 

.00004 

.00004 

.00005 

.00004 

.00004 

r 
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TABLE  B-6.  CRACK  OPENING  DISFIACEMENTS  FOR  A1  2024-T3  CENTER- CRACKED 
PANELS  FOR  A CYCLIC  LOADS  OBTAINED  FR(»1  AN  ELASTIC-PLASTIC 
FINITE  ELEMENT  ANALYSIS 


a > 0.52  inches  w ■ 12  Inches 

Crack  Opening  Displacement  (Inches)  At  Various 
Positions  On  Crack  Face  (Inches  from  Center^ 

x-.l  X-.2  X-.3  X-.4  X-.461 

CYCLE  2* 

Mean  Stress  ■ 13,125  psi  Stress  Amplitude  ■ 13,125  psl 


Loading  Phase 


3,281 

.00038 

.00037 

.00036 

.00033 

.00029 

.00027 

6,562 

.00069 

.00068 

.00065 

.00058 

.00048 

.00040 

9,844 

.00100 

.00099 

.00093 

.00083 

.00067 

.00052 

13,125 

.00131 

.00129 

.00122 

.00108 

.00086 

.00065 

16,406 

.00163 

.00160 

.00151 

.00133 

.00105 

.00077 

19,688 

.00195 

.00191 

.00179 

.00158 

.00124 

.00090 

22,969 

.00226 

.00222 

.00208 

.00183 

.00142 

.00102 

26,250 

.00257 

.00252 

.00237 

.00209 

.00161 

.00115 

22,496 

.00221 

Unloading  Phase 

.00217  .00204 

.00180 

.00139 

.00101 

18,746 

.00186 

.00182 

.00171 

.00151 

.00118 

.00086 

14,997 

.00150 

.00147 

.00138 

.00123 

.00097 

.00072 

11,248 

.00114 

.00112 

.00105 

.00094 

.00075 

.00058 

7,498 

.00078 

.00077 

.00073 

.00065 

.00054 

.00043 

3,749 

.00042 

.00042 

.00040 

.00037 

.00032 

.00029 

0 

.00006 

.00007 

.00007 

.00006 

.00011 

.00015 

* Results  for  Cycle  1 listed  in  Table  B. 


Applied 

Stress, 

psi 
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TABLE  B-7.  CKACK  OPEHING  DISPLACEMENTS  FOR  A1  2024-T3  CENTER- CRACKED 
PANELS  FOR  CYCLIC  STRESS  LOADING  OBTAINED  FRO!  AN  EIASTIC- 
PLASTIC  FINITE  ELEMENT  ANALYSIS  (Plane  Stress) 


a * 0.52  Inches  12  inches 

Crack  Opening  Displacement  (Inches)  At  Various 
Positions  On  Crack  Face  (Inches)  From  Center 

x-O  3P.1  X-.2  X-.3  X-.4  X-.461 


CYCLE  1 

Mean  Stress  - 19,688  psl  Stress  Amplitude  - 6,562  psi 


Loading  Phase 


22,875 

.00220 

.00215 

.00201 

.00176 

.00133 

.00090 

24,000 

.00232 

.00228 

.00213 

.00187 

.00142 

.00098 

25,125 

.00245 

.00240 

.00225 

.00198 

.00152 

.00170 

26,250 

.00257 

.00252 

.00237 

.00209 

.00161 

.00115 

21,872 

.00216 

Unloading  Phase 

.00211  .00198 

.00175 

.00136 

.00098 

17,498 

.00174 

.00170 

.00160 

.00142 

.00111 

.00082 

13,125 

.00132 

.00129 

.00122 

.00108 

.00086 

.00065 

16,406 

.00163 

.00160 

CYCLE  2 

Loading  Phase 

.00151 

.00133 

.00105 

.00077 

19,688 

.00195 

.00191 

.00179 

.00158 

.00124 

.00090 

22,968 

.00226 

.00222 

.00208 

.00184 

.00142 

.00102 

26,250 

.00264 

.00257 

.00237 

.00209 

.00161 

.00115 

Unloading  Phase 

21,872 

.00216 

.00211 

.00198 

.00175 

.00136 

.00098 

17,498 

.00174 

.00170 

.00160 

.00142 

.00111 

.00082 

13,125 

.00132 

.00129 

.00122 

.00108 

.00086 

.00065 

Applied 

Stress, 

psi 
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TABLE  B-8.  CBACK  OPENING  DISPLACEMENTS  FOR  A1  2024-T3  CENTER-CRACKED 
PANELS  FOR  CYCLIC  STRESS  LOADING  OBTAINED  FROM  AN  ELASTIC 
PLASTIC  FINITE  ELEMENT  ANALYSIS  (Plane  Stress) 


a * 0.52  inches  w ■ 12  Inches 

Applied  Crack  Opening  Displacements  (Inches)  At  Various 


Stress, 

psl 

Positions  On  Crack  Face  (Inches  From  Center) 

x-0 

x^.l 

x*.2 

x».3 

3P".4 

3^.461 

CYCLE  1 

Mean  Stress  ■ 

26,250  psl 

Stress  Amplitude  ••  13,125  psl 

Loadlns  Phase 

23,954 

.00232 

.00227 

.00212 

.00186 

.00142 

.00098 

26,156 

.00256 

.00251 

.00236 

.00208 

.00160 

.00114 

28,360 

.00282 

.00276 

.00260 

.00230 

.00180 

.00132 

30,563 

.00309 

.00303 

.00285 

.00254 

.00251 

.00152 

32,766 

.00337 

.00331 

.00313 

.00280 

.00225 

.00174 

34,968 

.00369 

.00362 

.00343 

.00309 

.00252 

.00198 

37,172 

.00403 

.00397 

.00377 

.00342 

.00283 

.00226 

39,375 

.00442 

.00435 

.00415 

.00379 

.00318 

.00258 

Unloadlns  Phase 

35,620 

.00404 

.00398 

.00381 

.00349 

.00296 

.00243 

31,871 

.00369 

.00364 

.00348 

.00320 

.00274 

.00229 

28,122 

.00333 

.00328 

.00315 

.00291 

.00253 

.00214 

24,373 

.00297 

.00293 

.00282 

.00263 

.00231 

.00200 

20,624 

.00261 

.00258 

.00249 

.00234 

.00210 

.00186 

16,874 

.00225 

.00223 

.00217 

.00205 

.00189 

.00171 

13,125 

.00190 

.00188 

.00184 

.00177 

.00167 

.00157 
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TABLE  B-8  (Continued)  CRACK  OPENING  DISPLACEMENTS  FOR  A1  2024-T3 

CENTER-CRACKED  PANELS  FOR  CYCLIC  STRESS  LOADING  OBTAINED  FROM 
AN  ELASTIC-PLASTIC  FINITE  ELEMENT  ANALYSIS  (Plane  Stress) 


Applied 
Stress , 
psi 

a ■ 

0.52  inches  w > 12  Inches 

Crack  Opening  DisplaceiDents  (Inches)  At  Various 
Positions  On  Crack  Face  (Inches  From  Center) 

x-0 

x*.l 

3^.2  x«.3 

x“.461 

CYCLE  2 

Mean 

Stress  ■ 

26,250  psi 

Stress  Amplitude  ■ 13,125  psi 

Loadins  Phase 

16,406 

.00220 

.00219 

.00212  .00202 

.00186 

.00169 

19,688 

.00252 

.00250 

.00241  .00227 

.00205 

.00182 

22,968 

.00284 

.00280 

.00270  .00252 

.00223 

.00195 

26,250 

.00215 

.00311 

.00299  .00277 

.00242 

.00207 

29,531 

.00346 

.00342 

.00327  .00302 

.00261 

.00220 

32,812 

.00378 

.00372 

.00356  .00327 

.00280 

.00232 

36,094 

.00409 

.00403 

.00385  .00352 

.00299 

.00245 

39,375 

.00441 

.00435 

.00414  .00378 

.00318 

.00257 

Unloadins  Phase 

35,620 

.60406 

,00400 

.00381  .00349 

.00296 

.00243 

31,871 

.00369 

.00364 

.00349  .00321 

.00275 

.00229 

28,122 

.00334 

.00329 

.00316  .00292 

.00254 

.00215 

24,373 

.00299 

.00294 

.00283  .00263 

.00232 

.00200 

20,624 

.00262 

.00259 

.00250  .00235 

.00211 

.00186 

16,874 

.00226 

.00224 

.00217  .00206 

.00189 

.00172 

13,125 

.00191 

.00189 

.00186  .00178 

.00170 

.00158 
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TABLE  B-8  (Continued^  CBACK  OPENING  DISPLACEMENTS  FOR  A1  2024-T3 

CENTER-CRACKED  PANELS  FOR  CYCLIC  STRESS  LOADING  OBTAINED  FROM 
AN  ELASTIC-PLASTIC  FINITE  ELEMENT  ANALYSIS  (Plane  Stress) 


a • 

0.52  Inches 

w ■ 12  Inches 

Applied 

Stress, 

Crack  Opening  Displacements  (Inches)  At  Various 
Positions  On  Crack  Face  (Inches  From  Center) 

psl 

3^0 

x-,1 

x*.2 

x*.3 

X-.4 

x“.461 

CYCLE  3 

Mean 

Stress  ■ 

26,250  psl 

Stress  Amplitude  - 13,125 

psl 

Loadlna  Phase 

16,406 

.00222 

.00220 

.00213 

.00202 

.00186 

.00167 

19,688 

.00253 

.00250 

.00242 

.00228 

.00225 

.00182 

22,968 

.00284 

.00281 

.00271 

.00253 

.00224 

.00195 

26,250 

.00316 

.00312 

.00209 

.00278 

.00243 

.00207 

29,531 

.00247 

.00342 

.00328 

.00303 

.00262 

.00220 

32,812 

.00378 

.00373 

.00357 

.00328 

.00280 

.00233 

36,094 

.00410 

.00404 

.00386 

.00353 

.00299 

.00245 

39,375 

.00441 

.00435 

.00414 

.00378 

.00318 

.00258 

Unloadlna  Phase 

35,620 

.00405 

.00399 

.00381 

.00349 

.00296 

.00243 

31,871 

.00369 

.00364 

.00349 

.00320 

.00275 

.00229 

28,122 

.00334 

.00329 

.00316 

.00292 

.00253 

.00215 

24,373 

.00298 

.00294 

.00283 

.00263 

.00232 

.00200 

20,624 

.00262 

.00259 

.00250 

.00235 

.00210 

.00186 

16,874 

.00226 

.00224 

.00217 

.00206 

.00189 

.00172 

13,125 

.00190 

.00189 

.00184 

.00177 

.00168 

.00157 

/ 
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TABLE  B*9.  CBACK  OPENING  DISPLACEMENTS  FOR  A1  2024-T3  CENTER- CRACKED 
PANELS  FOR  CYCLIC  STRESS  LOADING  OBTAINED  FROM  AN  EIASTIC- 
PLASTIC  FINITE  ELEMENT  ANALYSIS  (Plane  Stress) 


a * 0.52  Inches  w > 12  Inches 

Crack  Opening  Displacements  (Inches)  At  Various 
Positions  On  Crack  Face  (Inches  From  Center') 

xfO  x-.l  x*.2  X-.3  x“.4  XP.461 


Applied 
Stress s 
psi 


CYCLE  1 

Mean  Stress  ■ 26,250  psi  Stress  Amplitude  > 13,125  psi 


Loading  Phase 


23,954 

.00232 

.00227 

.00212 

.00186 

.00142 

.00098 

26,156 

.00256 

.00251 

.00236 

.00208 

.00160 

.00114 

28,360 

.00282 

.00276 

.00260 

.00230 

.00180 

.00132 

30,563 

.00309 

.00303 

.00285 

.00254 

.00201 

.00152 

32,766 

.00337 

.00331 

.00313 

.00280 

.00225 

.00174 

34,968 

.00369 

.00362 

.00343 

.00309 

.00252 

.00198 

37,172 

.00403 

.00397 

.00377 

.00342 

.00283 

.00226 

39,375 

.00442 

.00435 

.00415 

.00379 

.00318 

.00258 

Unloading  Phase 


35,620 

.00404 

.00398 

.00381 

.00349 

.00296 

.00243 

31,871 

.00369 

.00364 

.00348 

.00320 

.00274 

.00229 

28,122 

.00333 

.00328 

.00315 

.00291 

.00253 

.00214 

24,373 

.00297 

.00293 

.00282 

.00263 

.00231 

.00200 

20,624 

.00261 

.00258 

.00249 

.00234 

.00210 

.00186 

16,874 

.00225 

.00223 

.00217 

.00205 

.00189 

.00171 

13,125 

.00190 

.00188 

.00184 

.00177 

.00167 

.00157 

I 
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TABLE  B-9.  (Continued)  CRACK  OPENING  DISPLACEMENTS  FOR  A1  2024-T3 

CENTER- CRACKED  PANELS  FOR  CYCLIC  STRESS  LOADING  OBTAINED 
FROM  AN  ELASTIC-PLASTIC  FINITE  ELEMENT  ANALYSIS  (Plane  Stress) 


Applied 

Stress, 

psl 

a “ 

0.52  Inches  w « 12  Inches 

Crack  Opening  Displacements  (Inches)  At  Various 
Positions  On  Crack  Face  (Inches  From  Center) 

x*0 

X“.l 

xf.2 

x*.3 

X-.4 

x>.461 

CYCLE  2 

Mean  Stress 

=■  19,688  psl 

Stress  Aiiq>lltude  > 6,562,  psl 

Loadlns  Phase 

16,406 

.00221 

.00219 

.00212 

.00202 

.00186 

.00169 

19,688 

.00252 

.00250 

.00241 

.00227 

.00205 

.00182 

22,969 

.00284 

.00280 

.00270 

.00252 

.00223 

.00195 

26,250 

.00315 

.00311 

.00299 

.00277 

.00242 

.00207 

Unloadlns  Phase 

21,872 

.00273 

.00270 

.00260 

.00244 

.00217 

.00190 

17,498 

.00231 

.00229 

.00222 

.00210 

.00192 

.00174 

13,125 

.00190 

.00188 

.00184 

.00178 

.00167 

.00157 
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TABLE  B-10.  CRACK  OPENING  DISPLACEMENT  FOR  A1  2024-T3  CENTER- CRACKED 
PANELS  FOR  CYCLIC  STRESS  LOADING  OBTAINED  FRCM  AN  ELASTIC- 
PLASTIC  FINITE  ELEMENT  ANALYSIS  (Plane  Stress) 


a > 0.52  Inches  w ~ 12  Inches 

Crack  Opening  Displacements  (Inches)  At  Various 
Positions  On  Crack  Face  (Inches  From  Center) 

xfO  x*.1  3^.2  xf.3  x“.4  x*.461 

CYCLE  1 

Mean  Stress  ■ 19,688  psl  Stress  Amplitude  > 19,688  psi 


Applied 

Stress, 

psl 


Loading  Phase 


23,954 

.00232 

. .00227 

.00212 

.00186 

.00142 

.00098 

26,156 

.00256 

.00251 

.00236 

.00208 

.00160 

.00114 

28,360 

.00282 

.00276 

.00260 

.00230 

.00180 

.00132 

30,563 

.00309 

.00303 

.00285 

.00254 

.00201 

.00152 

32,766 

.00337 

.00331 

.00313 

.00280 

.00225 

.00174 

34,968 

.00369 

.00362 

.00343 

.00309 

.00252 

.00198 

37,172 

.00403 

.00397 

.00377 

.00342 

.00283 

.00226 

39,375 

.00442 

.00435 

.00415 

.00379 

.00318 

.00258 

Unloading  Phase 


35,620 

.00404 

.00398 

.00381 

.00349 

.00296 

.00243 

31,871 

.00369 

.00364 

.00348 

.00320 

.00274 

.00229 

28,122 

.00333 

.00328 

.00315 

.00291 

.00253 

.00214 

24,373 

.00297 

.00293 

.00282 

.00263 

.00231 

.00200 

20,624 

.00261 

.00258 

.00249 

.00234 

.00210 

.00186 

16,874 

.00225 

.00223 

.00217 

.00205 

.00189 

.00171 

13,125 

.00190 

.00188 

.00184 

.00177 

.00167 

.00157 

9,844 

.00158 

.00157 

.00135 

.00152 

.00148 

.00144 

6,562 

.00127 

.00127 

.00126 

.00127 

.00129 

.00132 

3,281 

.00095 

.00096 

.00098 

.00101 

.00111 

.00119 

0 

.00064 

.00065 

.00069 

.00076 

.00092 

.00107 
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TABLE  B-10.  (Continued)  CBACK  OPENING  DISPLACEMENT  FOR  A1  2024-T3 

CENTER-CRACKED  PANELS  FOR  CYCLIC  STRESS  LOADING  OBTAINED 
FROM  AN  ELASTIC-PLASTIC  FINITE  ELEMENT  ANALYSIS  (Plane  Stress) 


a ■ 

0.52  Inches 

w ■ 

12  Inches 

Applied 

Stress, 

Crack  Opening  Displacements  (Inches)  At  Various 
Positions  On  Crack  Face  (Inches  From  Center) 

psl 

x-0 

X“.l 

x«.2 

X-.3 

X-.4 

»«.461 

CYCLE  2 

Mean 

Stress 

* 19,688  psl 

Stress  Amplitude  « 19,688 

psl 

Loadlna  Phase 

3,281 

.00095 

.00096 

.00098 

.00101 

.00111 

.00119 

6,562 

.00127 

.00127 

.00126 

.00127 

.00129 

.00132 

9,844 

.00158 

.00157 

.00155 

.00152 

.00148 

.00144 

13,125 

.00190 

.00188 

.00184 

.00177 

.00167 

.00157 

16,406 

.00221 

.00219 

.00212 

.00202 

.00186 

.00169 

19,688 

.00252 

.00250 

.00241 

.00227 

.00205 

.00182 

22,969 

.00284 

.00280 

.00270 

.00252 

.00223 

.00195 

26,250 

.00215 

.00311 

.00299 

.00277 

.00242 

.00207 

29,531 

.00346 

.00342 

.00327 

.00302 

.00261 

.00220 

32,812 

.00378 

.00372 

.00356 

.00327 

.00280 

.00232 

36,094 

.00409 

.00403 

.00385 

.00352 

.00299 

.00245 

39,375 

.00441 

.00435 

.00414 

.00378 

.00318 

.00258 

Unloadlna  Phase 

35,796 

.00407 

.00401 

.00383 

.00351 

.00297 

.00244 

32,212 

.00373 

.00368 

.00352 

.00323 

.00277 

.00230 

28,633 

.00338 

.00334 

.00320 

.00296 

.00256 

.00217 

25,054 

.00304 

.00301 

.00289 

.00268 

.00236 

.00203 

21,475 

.00270 

.00267 

.00258 

.00241 

.00215 

.00189 

17,895 

.00236 

.00233 

.00226 

.00214 

.00195 

.00175 

U,316 

.00202 

.00200 

.00195 

.00186 

.00174 

.00162 

10,737 

.00167 

.00166 

.00164 

.00159 

.00154 

.00148 

7,158 

.00133 

.00133 

.00132 

.00132 

.00133 

.00134 

3,579 

.00099 

.00099 

.00101 

.00104 

.00113 

.00121 

0 

.00065 

.00066 

.00070 

.00077 

.00092 

.00107 
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TABLE  B-11.  CRACK  OPENING  DISPLACEMENTS  FOR  A1  2024-T3  CENTER- CRACKEt 
PANEL  FOR  A CYCLIC  STRESS  LQ4D  OBTAINED  FROM  AN  ELASTIC- 
PLASTIC  FINITE  ELEMENT  ANALYSIS  (Plane  Stress) 


a - 0.52  Inches  w - 12  Inches 

Crack  Opening  Displacements  (Inches)  At  Various 
Positions  On  Crack  Face  (Inches  From  Center) 

x»0  x-.l  X-.2  X-.3  x*.4  x«.461 


Applied 

Stress, 

psl 


CYCLE  1 

Mean  Stress  » 13,125  psl  Stress  Amplitude  > 13,125  psl 

Loading  Phase 


22,875 

.00220 

.00215 

.00201 

.00176 

.00133 

.00090 

24,000 

.00238 

.00232 

.00213 

.00187 

.00142 

.00098 

25,125 

.00245 

.00240 

.00225 

.00198 

.00152 

.00107 

26,250 

.00257 

.00252 

.00237 

.00209 

.00161 

.00115 

Unloading  Phase* 

22,500 

.00249 

.00245 

.00233 

.00212 

.00180 

.00151 

18,750 

.00208 

.00205 

.00196 

.00179 

.00152 

.00129 

15,000 

.00168 

.00166 

.00151 

.00145 

.00124 

.00106 

11,250 

.00128 

.00126 

.00121 

.00111 

.00096 

.00083 

7,500 

.00088 

.00087 

.00083 

.00077 

.00068 

.00061 

3,750 

.00048 

.00047 

.00046 

.00043 

.00040 

.00038 

0 

.00008 

.00008 

.00008 

.00008 

.00012 

.00015 

* Unloaded 

after  crack 

: extension 

of  0.06 

inches  at 

maximum  load. 

i t 

i » 
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TABLE  B-11.  (Continued)  CRACK  OPENING  DISPLACEMENTS  FOR  A1  2024-T3 
CENTER-CRACKED  PANEL  FOR  A CYCLIC  STRESS  LOAD  OBTAINED 
FROM  AN  ELASTIC-PLASTIC  FINITE  ELEMENT  ANALYSIS  (Plane  Stress) 


Applied 

Stress, 

psl 

a ■ 

0.52  Inches  w ■ 12  Inches 

Crack  Opening  Displacements  (Inches)  At  Various 
Positions  On  Crack  Face  (Inches  From  Center) 

x"0 

x«.l 

X-.2 

X-.3 

x“.4 

x“.461 

CYCLE  2 

Mean  Stress 

- 13,125  psl  Stress  Amplitude  > 13,  125  psl 

Loadlns  Phase 

3,281 

.00043 

.00043 

.00041 

.00039 

.00036 

.00035 

6,582 

.00078 

.00077 

.00074 

.00069 

.00061 

.00055 

9,844 

.00113 

.00112 

.00107 

.00098 

.00085 

.00075 

13,125 

.00148 

.00146 

.00139 

.00128 

.00110 

.00095 

16,406 

.00183 

.00181 

.00172 

.00157 

.00134 

.00114 

19,688 

.00218 

.00215 

.00205 

.00187 

.00159 

.00134 

22,969 

.00254 

.00250 

.00238 

.00217 

.00183 

.00154 

26,250 

.00289 

.00284 

.00271 

.00246 

.00208 

.00174 

Unloading  Phase* 

22,500 

.00277 

.00274 

.00263 

.00245 

.00218 

.00195 

18,750 

.00233 

.00230 

.00221 

.00207 

.00184 

.00166 

15,000 

.00188 

.00186 

.00179 

.00168 

.00150 

.00137 

11,250 

.00144 

.00142 

.00137 

.00129 

.00117 

.00107 

7,500 

.00100 

.00099 

.00095 

.00090 

.00083 

.00077 

3,750 

.00055 

.00055 

.00053 

.00051 

.00049 

.00048 

0 

.00011 

.00011 

.00011 

.00012 

.00015 

.00018 

* Unloaded  after  crack  extension  of  0.06  Inches  at  bmxIioux  load 
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TABLE  B-11.  (Continued)  CRACK  OPENING  DISPLACEMENTS  FOR  A1  2024-T3 
CENTER-CRACKED  PANEL  FOR  A CYCLIC  STRESS  LOAD  OBTAINED 
FRW  AN  ELASTIC-PLASTIC  FINITE -ELEMENT  ANALYSIS  (Plane  Stress) 


Applied 

Stress, 

a ■ 

0.52  Inches  w > 12  Inches 

Crack  Opening  Displacements  (Inches)  At  Various 
Positions  On  Crack  Face  (Inches  From  Center) 

psl 

x«0 

XF.l 

X-.2 

x»,3 

X-.4 

x«.461 

Mean 

Stress 

= 13,125  psl 

CYCLE  3 

Stress  Amplitude  > 13, 

125  psl 

3,281 

.00050 

Loading  Phase 

.00049  .00048 

.00046 

.00044 

.00044 

6,562 

.00089 

.00087 

.00085 

.00080 

.00074 

.00070 

9,844 

.00127 

.00126 

.00122 

.00114 

.00104 

.00096 

13,125 

.00166 

.00164 

.00158 

.00148 

.00133 

.00122 

16,406 

. 00205 

.00203 

.00195 

.00182 

.00163 

.00148 

19,688 

.00244 

.00241 

.00232 

.00216 

.00193 

.00174 

22,969 

.00283 

.00279 

.00269 

.00250 

.00222 

.00199 

26,250 

.00322 

.00318 

.00305 

.00284 

.00252 

.00225 

22,500 

.00304 

Unloading  Phase 

.00300  .00291 

k 

.00275 

.00251 

.00233 

18,750 

.00255 

.00252 

.0C245 

.00231 

.00212 

.00197 

15,000 

.00207 

.00204 

.00198 

.00188 

.00173 

.00161 

11,250 

.00158 

.00156 

.00152 

.00144 

.00133 

.00126 

7,500 

.00109 

.00108 

.00105 

.00100 

.00094 

.00090 

3,750 

.000604 

. 000600 

.000588 

.000569 

.000549 

.000542 

0 

.000117 

. 000118 

. 000123 

.000132 

.000156 

.000185 

Unloaded  after  crack  extension  of  0.06  Inches  at  maximum  load. 
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TABLE  B-11.  (Continued)  CRACK  OPENING  DISPIACEMENTS  FOR  A1  2024-T3 
CENTER- CRACKED  PANEL  FOR  A CYCLIC  STRESS  LOAD  OBTAINED 
FRO!  AN  ELASTIC-PLASTIC  FINITE- ELEMENT  ANALYSIS  (Plane  Stress) 


a *= 

0.52  Inches 

W “ 

12  inches 

Applied 

Stress, 

Crack  Opening  Displacements  (Inches)  At  Various 
Positions  On  Crack  Face  (Inches  From  Center) 

psl 

x=0 

x-.l 

x*.2 

JC“.3 

IP*. 4 

x=.461 

CYCLE  4 

Mean 

Stress 

= 13,125  psl 

Stress  Amplitude  ^13, 

125  psl 

Loadlne  Phase 

3,281 

.00054 

.00053 

.00053 

.00051 

.00050 

.00049 

6,562 

.00097 

.00096 

.00094 

.00090 

.00084 

.00081 

9,844 

.00140 

.00138 

.00134 

.00128 

.00119 

.00112 

13,125 

.00182 

.00180 

.00175 

.00166 

.00153 

.00143 

16,406 

.00225 

.00222 

.00216 

.00204 

.00187 

.00175 

19,688 

.00267 

.00264 

.00256 

.00242 

.00222 

.00206 

22,969 

.00310 

.00307 

.00297 

.00281 

.00256 

.00237 

26,250 

.00353 

.00349 

.00338 

.00320 

.00291 

.00268 

* 

Unloadlns  Phase 

22,500 

.00332 

.00329 

.00321 

.00306 

.00285 

.00269 

18,750 

.00279 

.00277 

.00270 

.00258 

.00241 

.00228 

15,000 

.00226 

.00224 

.00219 

.00209 

.00196 

.00186 

11,250 

.00173 

.00172 

.00168 

.00161 

.00151 

.00145 

7,500 

.00120 

.00119 

.00117 

.00112 

.00107 

.00103 

3,750 

.000670 

.000666 

.000656 

.000639 

.000623 

.000620 

0 

.000139 

.000141 

.000145 

.000155 

.000178 

.000206 

* " - . — 

Unloaded  after  crack  extension  of  0.06  Inches  at  load. 
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DISTRIBUTION  LIST 


REPORT  NO.  NADC-77255-30 
AIRTASK  NO.  WR  02303001,  WORK  UNIT  DG-602 

No.  of 

NAVAIRSYSCOM,  AIR -50174  Co£ie8 

(2  for  retention,  2 for  AIR-530,  1 for  AIR-530215,  1 for  AIR-530221C, 

2 for  A1R-320B) 8 

NAVAIRTESTCEN,  Patuxent  River,  Maryland  1 

NAVAVNSAFECEN,  NAS,  Norfolk,  Virginia  1 

CNAVANTRA,  NAS,  Corpus  Chris ti,  Texas  1 

CNABATRA,  NAS,  Pensacola,  Florida  1 

CNARESTRA,  NAS,  Glenview,  Illinois  1 

CNATRA,  NAS,  Pensacola,  Florida  1 

NAVAIRSYSCOMREPLANT  1 

NAVAIRSYSCOMREPCENT  1 

NAVAIRSYSCOMREPAC  1 

NAVAIREWORKFAC,  NAS,  Alameda,  California  1 

NAVAIREWORKFAC,  NAS,  Jacksonville,  Florida  1 

NAVAIREWORKFAC,  NAS,  Norfolk,  Virginia  1 

NAVAIREWORKFAC,  NAS,  Pensacola,  Florida  1 

NAVAIREWORKFAC,  NAS,  Quonset  Point,  Rhode  Island  1 

NAVAIREWORKFAC,  NAS,  San  Diego,  California  1 

NAVAIREWORKFAC,  NAS,  Cherry  Point,  North  Carolina  1 

COMNAVAIRLANT 1 

COMNAVAIRFAC  1 

NWL,  Dahlgren,  Virginia  (Attention  Mr.  Morton)  ...  1 

USAF  Systems  Command,  WPAFB,  Ohio  45433 

Attention  FBR 1 

Attention  FB 1 

Attention  LLD 1 

Attention  SEFS I 

Attention  FYA 1 

Attention  LAM 1 

Attention  FBA  1 

Attention  LPH  ............  I 

USA  AMMRC,  Watertown,  Massachusetts  1 

USA  APG,  Aberdeen,  Maryland  I 

USA  AMRDL,  Fort  Eustis,  Virginia  (Attention  Mr.  Berrisford)  . 1 

USA  AVSCOM,  St.  Louis,  Missouri  (Attention  AMSAV-GR)  1 

Defense  Research  and  Development  Staff,  British  Embassy, 

Washington,  D.C.,  via  NAVAIR  (AIR-5302)  I 

Canadian  Joint  Staff,  Navy  Member, 

Washington,  D.C.,  via  NAVAIR  (AIR-5302)  I 

Technical  Advisory,  AFIAS-B,  Directorate  of  Aerospace  Safety, 

Morton  AFB,  California I 

DDC 12 

NAVSEASYSCOM,  Washington,  D.C.  20362  (Attention  Mr,  C,  Pohler, 

Code  035) • 1 

NAVSHIPRADCEN,  Bethesda,  Maryland  20034 

(Attention  Mr.  A.  B.  Stavovy  730) 1 

Bell  Acrosystem  Co.,  Buffalo,  New  York  14205  1 

Bell  Helicopter  Co.,  Fort  Worth,  Texas  76101  I 

Boeing  Co.,  Airplane  Div.,  SeattU,  Washington  98124 

(Attention  Mr.  T.  Porter)  1 
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Boeing  Co.,  Airplane  Div. , Wichita,  Kansas  67210  1 

Boeing  Co.,  Vertol  Div.,  Philadelphia,  Pennsylvania  19142  1 

McDonnell  Douglas  Aircraft  Corp.,  Aircraft  Div., 

Long  Beach,  California  90801  1 

General  Dynamlcs/Convalr,  San  Diego,  California  92112  1 

General  Dynamics  Corp.,  Fort  Worth,  Texas  76101  1 

Goodyear  Aerospace  Corp.,  Akron,  Ohio  44305  1 

Grumman  Aerospace  Corp.,  Bethpage,  Long  Island,  New  York  11714  ....  1 

Aircraft-Missiles  Div.,  Falrchlld-Hlller  Corp 1 

Hagerstown,  Maryland  21740  ..  1 

Kaman  Aircraft  Corp.,  Bloomfield,  Connecticut  06002  1 

Lockheed  Aircraft  Corp.,  Lockheed-Callfornla  Co., 

Burbank,  California  91503  1 

Lockheed  Aircraft  Corp.,  Lockheed -Georgia  Co., 

Marietta,  Georgia  30061  1 

LTV  Aerospace  Corp.,  Dallas,  Texas  75222  1 

Martin  Co.,  Baltimore,  Maryland  21203  1 

McDonnell  Douglas  Aircraft  Corp.,  St.  Louis,  Missouri  63166  1 

Rockwell  International,  Columbus  Aircraft  Div., 

Columbus,  Ot^lo  43216  (Attention  Mr.  0.  Acker) 1 

Rockwell  International,  Los  Angeles,  California  90053 

(Attention  Mr.  G.  Fitch) 1 

Northrop  Corp.,  Aircraft  Div.,  Hawthorne,  California  90250  1 

Republic  Aviation  Div.,  Falrchlld-Hlller  Corp., 

Farmlngdale,  Long  Island,  Mew  York  11735  1 

Sikorsky  Aircraft  Co.,  Stratford,  Connecticut  06497  . . 1 

Hdqtrs.,  R&T  Div.,  AFSC,  Bolling  AFB,  Washington,  D.C.  22209  1 

Office  of  Aerospace  Research,  Arlington,  Virginia  22209  1 

FAA  (FS-120),  Washington,  D.C.,  20553  1 

Scientific  and  Technical  Information  Facility, 

College  Park,  Maryland  20740  (NASA  Rep. ) 2 

Administrator,  NA^,  Washington,  D.C.  20546  1 

NASA-Langley  Research  Center,  Materials  Div.,  HaBq>ton  Virginia  23365 

(Attention  Mr.  H.  F.  Hardrath) 1 

National  BuStds,  Washington,  D.C.  20234  1 

Office  of  Naval  Research,  Washington,  D.C.  20362 

(Attention  Dr.  N.  Perrons) 1 

Director,  Naval  Research  Lab,  Washington,  D.C.  20390  1 

Midwest  Research  Institute,  lUnsas  City,  Missouri  64110  1 

University  of  Illinois,  Urbana,  Illinois  61803  (Attention  Prof.  T.  J. 

Dolan  and  Prof.  J.  Morrow)  I each 

University  of  Kansas,  Lawrence,  Kansas  66044  1 

University  of  Michigan,  Ann  Arbor,  Michigan  48105  .....  1 

University  of  Minnesota,  Minneapolis,  Minnesota  55455  1 

Alcoa,  Alcoa  Labs,  Alcoa  Center,  Pennsylvania  15069 

(Attention  Mr.  J.  G.  Kaufman) 1 

Lehigh  University,  Bethlehem,  Pennsylvania  18015 

(Attention  Prof.  6.  C.  Sih) 1 

Balfoor-Stulan,  Inc.,  Traverse  City,  Michigan  49684  1 

Cornell  Aero.  Imb,  Buffalo,  Mew  York  14221  1 
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I Metals  & Ceramics  Information  Center,  Battelle,  Columbus  Laboratories 

^ Columbus,  Ohio  43201  1 

NASA,  Lewis  Research  Center,  Cleveland,  Ohio  44153 

(Attention  Technical  Library)  1 

Naval  Postgraduate  School,  Monterey,  California  93940 

(Attention  Prof.  Lindsey)  1 

Battelle,  Pacific  Northwest  Laboratories,  P.  0.  Box  999,  Richland, 

Washington  99352  (Attention  Mr.  W.  E.  Anderson)  1 

School  of  Engineering,  George  Washington  University,  Washington, 

D.C.  20006  (Attention  Prof.  H.  Liebowitz).  1 

Dept,  of  Metallurgy,  University  of  Connecticut,  Storrs,  Connecticut 

06268  (Attention  Prof.  A.  J.  McEvlly)  1 

Dept,  of  Mech.  Eng.,  Washington  University,  St.  Louis,  Missouri 

(Attention  Prof.  P.  C.  Paris)  1 

Materials  Dept.,  6531  Boelter  Hall,  University  of  California  at  Los 
Angeles,  Los  Angeles,  California  90024  (Attention  Prof. 

A.  S.  Tetelman) 1 

National  Transportation  Safety  Board,  800  Independence,  Washington, 

D.C.  20594  (Attention  Mr.  W.  D.  Cowan)  1 

Pratt  and  Whitney  Aircraft,  P.  0.  Box  611,  Middletown,  Connecticut 

* (Attention  Mr.  C.  A.  Rau) 1 

Aeronautics  Dept.,  Delft  University,  Delft,  Holland  (Attention 

Prof.  J.  Schljve) 1 

' FAA,  800  Independence,  Washington,  D.  C.  20594  (Attention  Mr.  D.  Kemp)  . 1 

Air  Force  Flight  Dynamics  Lab,  WPAFB,  Ohio  45433  (Attention 

Mr.  R.  M.  Engle) 1 

VanderBllt  University,  F.  0.  Box  3245  Station  B,  Nashville, 

Tennessee  37235  (Attention  Prof.  P.  Packman) . 1 

Air  Force  Flight  Dynamics  Lab.,  Wright-Patterson  Air  Force  Base, 

Ohio  45433,  (Attention  Mr.  H.  A.  Wood) 1 

I NASA-Langley  Research  Center,  Hampton,  Virginia  23665 

'•  (Attention  Dr.  J.  C.  Newman,  Jr.) 1 

Max-Planck-Instltut  for  Elsenforschung,  4000  Dusseldorf, 

I Max-Planck-Strasse  1,  Germany  (Attention  Prof.  P.  Neumann)  1 

i Lehigh  University,  Metallurgy  and  Materials  Science,  Bethlehem, 

Pennsylvania  18015  (Attention  Prof.  R.  W.  Hertzberg) 1 

Southwest  Research  Institute,  8500  Culebra  Rd.,  San  Antonio, 

Texas  78228  (Attention  Dr.  J.  Lankford)  1 

Electric  Power  Research  Institute,  P.O.  Box  10412,  3412  Hlllvlew 

Avenue,  Palo  Alto,  California  94303  (Attention  R.  Stahlkopf).  ...  1 

The  University  of  Newcastle  Upon  Tyne,  Newcastle  Upon  Tyne  NEl  7RU 

England  (Attention  Dr.  J.  T.  Evans,  Dept,  of  Metallurgy)  1 

United  Kingdom  Atoaiic  Energy  Authority,  11  Charles  2 St.,  SWl, 

. London,  England  (Attention  B.  Tomkins)  . • 1 

Purdue  University,  Lafayette,  Indiana  47907  (Attention  B.  M.  Hillberry).  1 
Massachusetts  Institute  of  Technology,  77  Massachusetts  Ave., 

Ca^rldge,  Massachusetts  02139  (Attention  R.  M.  Pelloux) 1 

U.S.  Naval  Research  Laboratory,  4555  Overlook  Ave.,  SW,  Washington,  D.C. , 

20032,  (Attention  T.  W.  Crooker)  1 
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University  of  Harvard,  Cambridge,  Massachusetts  02138 

(Attention  J.  W.  Hutchinson) 1 

U.S.  Nuclear  Regulatory  Commission,  Washington,  D.C.,  20555 

(Attention  R.  E.  Johnson) 1 


University  of  Iowa,  Iowa  City,  Iowa,  52252  (Attenslon  R.  I.  Stephens).  . ^ \ 


» 


